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I. INTRODUCTION 


The design, development and flight tests of the Digital Integrated Auto- 
matic Landing System (DIALS) described in this report were performed within 
the context of the Advanced Transport Operating Systems (AT OPS) program con- 
ducted by the Langley Research Center of the National Aeronautics and Space 
Administration. ATOPS is a research and development program aimed at devel- 
oping capabilities for increased terminal area capacity, safe and accurate 
flight in adverse weather conditions including shear winds, noise reduction, 
the avoidance of wake vortices and reduced fuel consumption. Advances in 
digital flight computers and modern control theory, coupled with accurate 
guidance information such as that provided by the Microwave Landing System 
can be effectively used to achieve some of these goals. The work described 
in this report has evolved within the ATOPS program framework [1]. W. 

The design of the Digital Integrated Automatic Landing System (DIALS) 
was completed in 1980. The automatic landing system was implemented on the 
Transport System Research Vehicle (TSRV) , a Boeing 737-100. DIALS was flight 
tested with successful test results by NASA’s Langley Research Center at the 
NASA Wallops Flight Center, demonstrating the application of modern control 
theory to a complex aircraft control design problem. The DIALS flight tests 
were completed in December 1981. 

Exceptionally few automatic control systems designed with modern optimal 
control techniques have been successfully flight tested. In fact, the author 
is aware of no such flight tests on jet aircraft prior to, or since, the 
DIALS flight tests. DIALS is also the first digital automatic landing system 
designed with a modern control structure and methodology which has been suc- 
cessfully flight tested, to the author’s knowledge. Thus, DIALS represents 
a milestone in the actual demonstration that linear modern control theory. 



when appropriately interpreted and applied, can result in working designs for 
the most complex systems and tasks. This link between theory and practice is 
one of the basic underlying objectives of the work performed. This report de- 
scribes the design of the Digital Integrated Automatic Landing System, and 
presents an analysis of the flight test results. 

During the DIALS flight tests, ten "hands-off” automatic landings were 
performed on the TSRV aircraft. These automatic landings were performed under 
diverse wind conditions, including a wind shear case. The final approach 
flights included the modes of localizer capture, steep glideslope capture, 
localizer and glideslope track, crab/decrab, and flare to touchdown. On 
three of the approaches, a conventional 3° glideslope was selected. On the 
remaining seven automatic landings, steep approaches with a 4.5° glideslope 
were made. Also, several approaches capturing and tracking 5° glideslopes 
were made without continuing till touchdown. 

DIALS uses the Microwave Landing System (MLS) which is a guidance system 
providing high accuracy position information in the form of azimuth, elevation, 
and range measurements. As the Microwave Landing System is less sensitive to 
weather conditions than conventional systems, automatic landing systems using 
the MLS can be used to reduce the congestion in terminal areas due to adverse 
weather conditions. Furthermore, the volumetric coverage provided by the MLS 
enables the use of curved flight paths and steep glideslopes in the final ap- 
proach and landing phase of the flight. The MLS guidance system consists of a 
DME providing range information, an azimuth antenna, generally colocated with 
the DME antenna, providing the aircraft Vs azimuth angle relative to the run- 
way centerline up to ±60°, and an elevation antenna located at the glidepath 
intercept point but offset to the side of the runway providing the aircraft’s 
elevation angle up to 20°. An onboard MLS receiver provides high accuracy 
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position information that can be used for steep approaches and curved flight 
paths in the terminal area* 

The control system uses MLS position information, as well as onboard 
sensor measurements but does not require inertial platforms. The system was 
developed using modern digital control methodologies [3] - [ 7] . The phases 
of the final approach and landing considered are localizer and glideslope 
capture and track, crab/decrab, and flare. The system captures, tracks, and 
flares from a steep glideslope selected by the pilot prior to engaging the 
capture mode. Thus, the control system modes designed are: 

1) Localizer capture 

2) Steep glideslope capture 

3) Localizer track 

4) Steep glideslope track 

5) De.crab 

6) Flare 

Two features of this sytem are: (1) the simultaneous capture of the 

localizer and glideslope and (2) the selectable steep glideslope (2.5° - 5.5° 
for the B-737). 

Typically, the aircraft is positioned, manually or automatically, below 
the selected glideslope so as to intercept the localizer at a selected air- 
speed, as shown in Figure 1. When the localizer or glideslope capture criteria 
are satisfied, the corresponding capture mode is engaged; so that depending on 
the initial aircraft position and attitude, the localizer and glideslope aan 
be captured simultaneously or sequentially. As the aircraft reaches the local- 
izer or glideslope, the corresponding track (or hold) modes are engaged. When 
a cross-wind component is present, the control system crabs the aircraft into 
the wind, followed by the decrab maneuver when the decrab altitude is reached. 
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A flare path which depends on the glideslope selected for a given approach 
is generated on-line, and the aircraft is controlled about this path until 
touchdown. 

The main considerations in the development of DIALS include: 

• low overshoots of the localizer and steep glideslope under 
adverse wind conditions 

• quick settling on the trajectory 

• overall smoothness of the flight during the final approach 

• accurate tracking in adverse weather conditions including 
gusts and wind shear. 

The ability to follow various steep glideslopes (selectable until glide- 
slope capture) provides a flexibility which can increase the efficiency of 
terminal area operations, reduce the noise perceived on the ground, and can 
be used for vortex avoidance when following a large aircraft, while reducing 
fuel consumption during a steep final approach. Capturing the glideslope 
and localizer simultaneously and with quick settling times allows close-in 
captures, while low overshoots of the localizer under adverse wind conditions 
enhance the independence of close parallel runway operations. Finally, the 
low degradation of the MLS information accuracy in adverse weather conditions 
enhances performance capabilities under low visibility conditions. 

The overall objective in the development of DIALS has been the direct- 
digital-design, using modern control methods, of an automatic landing system 
which controls the aircraft for glideslope angles ranging between 2.5° - 5.5°, 
nominal airspeeds between 115 - 135 knots, aircraft weights between 70,000 - 
90,000 lb, c.g., locations between .2 — .3 under adverse wind conditions and 
low visibility. As the control objectives during different phases of the 
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final approach vary, the most important desirable characteristics of the 
control law in the various portions of the flight are different. For example, 
while low overshoot characteristics are important during capture, vertical 
path accuracy becomes more significant in flare. A more detailed description 
of the desirable characteristics in the various control modes is given in 
Section IV. To achieve these varying objectives, it appears that some changes 
in the control law, from one phase of the flight to the next, are necessary. 
From a strict optimal control point of view, it would be necessary to change 
all the control gains for each phase of flight and aircraft condition. To 
avoid many changes, the following approach was adopted: 

1) obtain an optimal control law for a nominal condition as 
a starting point, 

2) using realistic nonlinear simulations, modify and update 
this law to account for the various nonlinearities in the 
actuator servomechanisms, hydraulic systems, engine dynam- 
ics and the aerodynamic response, as well as to enhance 
desirable characteristics which may not adequately be re- 
flected in a quadratic cost function, 

3) identify a small number of modifications which can be made 
at each phase of flight to achieve the differing objectives 
of each mode. 

The employment of modern digital design techniques is well-suited to the 
discrete nature of the MLS position information and the use of digital flight 
computers. From the point of view of using lower sampling rates, it is also 
preferable to the alternate approach of designing an analog system whose re- 
sponse is then approximated by a digital system, as the latter approach usually 
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results in rather high sampling rates. The control law was developed using 


optimal stochastic digital control methods. 

Section II summarizes the design of the control law. Section III de- 
scribes the implementation of DIALS on the TSRV flight computer. Section IV 
describes the modes of the automatic landing system and presents the flight 
test results. 
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II* CONTROL LAW FORMULATION 


The design of the Digital Integrated Automatic Landing System (DIALS) is 
based on modern optimal control design techniques. The design uses a sampled- 
data formulation of the Linear-Quadrat ic-Gaussian (LQG) [ 8 ] problem coupled 
with optimal disturbance accommodation [7] . The basic structure of the 
control system contains a fixed gain (steady-state) Kalman filter which esti- 
mates the aircraft states as well as the wind velocities, a desired flight 
path generator which produces the flight path* and attitude commands as well 
as the offsets or errors from the commanded trajectory during the various 
phases of the final approach and landing, and the control command unit con- 
taining the control gains, mode switching, limiters, etc. 

The approach is to develop a linear design model for the plant consist- 
ing of the aircraft dynamics, wind conditions, actuator dynamics, and the 
control law structure desired. Using the design model, a digital control law 
to track the desired trajectory is obtained by applying optimal stochastic 
sampled-date techniques. The control law obtained is tested in a detailed non- 
linear simulation of the aircraft and its subsystems. Finally, the control 
law is flight tested in automatic landings to test the design method, with the 
performance measure of the design method being the capability of the aircraft 
to perform the various functions accurately and completely. 

During most of the final approach and landing phases of flight, the air- 
craft maintains a level wings attitude with small deviations of bank angle. 

Thus, the flight condition corresponding to the tracking of the desired glide- 
slope and the localizer was selected as the reference condition about which the 
perturbation equations were obtained. In this condition, the longitudinal and 
lateral perturbation equations are decoupled; i.e. , small changes in the lateral 
variables affect the longitudinal variables only as second order effects; 
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conversely, the effect of the longitudinal variables on the lateral motion is 
also of second order. Thus, the control and modeling of the longitudinal and 
lateral dynamics are considered separately. 

Three sets of coordinate axes are used: the Earth-fixed axes, the body 

axes, and the stability axes. The Earth-fixed coordinate frame (x^, y , z^) 
has its origin fixed at the glidepath intercept point (GPIP) on the runway. 

The x & axis is along the runway centerline, the direction in which the aircraft 

* 

lands being chosen positive along the x^. The axis is the local vertical, 
positive downwards; y^ is perpendicular to both x^ and z^, with its positive 
end directed so as to make the coordinate frame right-handed. The Earth is 
assumed to be stationary with respect to inertial space; so that the Earth- 
fixed axes form an inertial frame. 

The body axes (x^, y^, z^) and the stability axes (x g , y g , z g ) are fixed 

to the body of the aircraft; i.e. , they are body-fixed axes. The origin of 

both axes is fixed at the aircraft center of mass. The x, axis is along the 

b 

fuselage reference line of the aircraft, positive towards the nose, the y^ axis 

is positive towards the tip of the right wing, the z^ axis is perpendicular to 

both x^ and y^ and is positive downwards (when the aircraft pitch angle is zero) . 

The stability axes (x , y , z ) are obtained from the body axes by a rotation 

s s s 

of a^, the steady-state angle of attack, about the y^ axis. The three sets of 
coordinate frames are shown in Figure 2. 

A. LONGITUDINAL DESIGN MODEL 

The general equations of motion for rigid aircraft can be linearized about 
a steady flight condition as described in (jl] , [lo] . The nominal flight condition 
used here corresponds to flight among the selected glideslope and localizer, at 
a constant airspeed with the flaps at 40 degrees and the gear down in the landing 
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conf iguration* The longitudinal perturbation equations for describing the air- 
craft’s motion in the vertical plane can be expressed in the stability axes as: 


mu = - mg cosy o 0 + q Q S (-C^, + 2C D + C T + 20^ ) u' 

( O X xo 


+ (C_ - CL ) a - G_ « 6e - CL* 6 s 4- c Xrr 6T 

Lo Da — D6e Dds T 6T 

— X 


m(w - U Q q) = - mg siny Q 9 + q Q S j - (C^, + 20^) u' 


- (C_ + C ) a - C T . a - C T q 

La Do — La — L£ • a - 


I yy q q o Sc | ^Mu’ + 2< W + ^ C Ma + C MTa^ - + C Ma - 

+ %^ +C MSe 4e + C M6s 4s + C M6T 61 j' 

u + u 
, w 

± = “u~ 

o 


a = a + a , 

— w 


SL = 


q + q. 


W 
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( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


where U is the reference inertial speed in the x direction, y the desired 
o s o 

glideslope angle, m the aircraft mass, 0 the perturbation in pitch angle, q the 
pitch rate, u the perturbation in inertial speed along the x g direction, w the 
inertial speed along the z g direction; u^, and q^ are the components due 
to wind, q Q is the steady-state value of the dynamic pressure at the selected 
airspeed, S is the effective wing area, 6e, 6s, and 6T are the perturbations of 
elevator, stabilizer and thrust, respectively. 

Equations (1) - (3) describe the linear and angular velocities of the air- 
craft in the vertical plane. The position of the aircraft can be obtained by 
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integrating the inertial velocity components over time. Thus, 


x; = ^ = L es (1,1) (1 + u') + l es (1,2)3 + l es (l,3)a , (7) 

O 

Z 

z; = I s = L es (3,1) (1 + u’) + L (3,2)3 + L es (3,3)a , (8) 

o 

where L (i, j ) is the element in i^ row and column of the matrix, L^g, 

which represents the transformation from the stability axes to the Earth-fixed 
axes. Note that a and 3 are inertial quantities and correspond to normalized 
velocity components in the stability axes; under no wind conditions these would 
be the same as the aerodynamic angle of attack and sideslip. The position 
equations can be rewritten in the following form 

x f = - siny 0 + cosy u 1, ■+ siny a + r] , (9) 

e r o Vo * o x 9 

= - cosy Q 0 - siny Q u f + cosy Q a + , (10) 

where 


n x = L es (1,1) (1 + u') - cosy o u' + siny o 9 + L Eg (1,2)3 

+ (L Eg (1,3) - slny o )a , (11) 

n z = Lg S (3,1) (1 + u') + siny Q u' + cosy o 0 + L Eg (3,2)3 

+ (L es (3,3) - cosY o )a , (12) 

In this form, the differential equations are linear with a forcing func- 
tion that contains the nonlinear part, which are second order terms with 
respect to the steady flight conditions considered. With this approach it is 
possible to use linear theory in the development of the filter and control law 
without neglecting the nonlinear terms completely. To obtain the second order 
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terms for the remaining variables, consider the equations 


0 = cos(}) q - sine}) r , 

(13) 

a = u + qw - rv , 

s 

(14) 

a ~ w - q(U 4- u) + pv 

z n o r 

s 

(15) 

Rearranging these equations, we obtain 


0 = q + P 0 , 

(16) 

a 

u ’ = u + V ’ 

o 

(17) 

a 

. z 

w s , , _ 

a = u = u + q + n a , 

(18) 


o o 

where 

ru = (coscj) - 1) q - sincj) r , (19) 

C7 


p , = - q a + r (3 , (20) 

u 

q = q u’ - p $ (21) 

a H 

It should be noted that the perturbation equations (1) - (3) neglect the 
second order terms represented by Pq, p^, , and p^. The second order terms 
introduce some coupling between the longitudinal and lateral equations. 

A simple first order model is used to represent the engine dynamics which 
are relatively slow due to spool-up time. The engine dynamics tend to behave 
with a 1 sec time constant when reducing thrust, while a 2 sec time constant 
behavior seems to dominate when increasing thrust, except when a rate limit 
is active to avoid enging overheating. The design model uses the 2 sec time 
constant to avoid lags in the actual operation of the control law. 
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ST = .5 ST + .298 Sth 


( 22 ) 


where ST is the thrust perturbation in units of 1000 lb per unit ST, and Sth 
is the throttle perturbation in throttle degrees. 

A control structure of rate command is used for the throttle and stabilizer 
to obtain relatively slow control commands. For this purpose, the design model 
is augmented by the equations. 


Sth = u 2 , 

(23) 

6s = , 

(24) 


where 5s is the stabilizer perturbation in radians. As the elevator actuator 
system is fast, with an approximate time constant of 1/16 sec, the elevator 
dynamics are neglected in the design model. 

To complete the model of the longitudinal equations of motion, it is 
necessary to model the wind velocities which affect the motion of the aircraft. 

The longitudinal wind model contains the components of steady wind velocities, 
turbulence and shear winds in the longitudinal axes. The turbulence model uses 
the Dry den spectra for the various components varying with altitude. The turbu- 
lence model has three components: u^ in the x^ direction, in the z^ direction, 

and q^ which models the effect of turbulence on the pitch rate of the aircraft. 

The u 1 component is independent of a and q ; however, a and q are correlated. 

8 § § 8 8 

The Dryden spectra can be factored using spectral factorization methods to obtain 
a linear system driven by white noise which generates an output having the desired 
spectral characteristics [9], [6], [ll]. Thus, the following transfer functions 
(i.e. , the ratio of the Laplace transforms of the gust generating system output 


to its input) are obtained to generate u f , a and q . 

§ 8 8 


G u (s) 


l +c Hs 

o 
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G (s) = 
a 


1 + 3 s 
o 


l + 2^s + 
o 


L \2 0 

w \ 2 

k) s 


( 26 ) 


q (s) s 

G ( s ) = = - 

q « 1 + 4r- s 


a (s) 
g 


(27) 


7TU 


where (x^ is the input to the system G (s) which generates q with the specified 

q g 

spectrum and cross-spectral density, and L , L are the scales of turbulence. 

a w 

The steady and shear wind in the longitudinal directions is modeled by 


sh 


sh 


’£3 


(28) 


% * 5 4 


(29) 


Thus, to simulate a specified shear profile for u\ with appropriate 

s 

initial conditions, e.g., to obtain a linear profile u T changing at a rate of 

s 

u\ , the initial condition for u T , is set to u f , and £ no is set equal to 
sho 5 sh sho £3 

zero; alternately, an impluse in will also achieve the same profile. The 
aircraft’s longitudinal equations of motion, the position equations, the actuator 
equations and the wind model above can be combined and after some manipulation 
can be expressed in state variable form. 


x-g = A q*q + B-u- + DfWjy + 

A/ A/ A/ Af A / A/ A/ A/ 


”2 A wl W 2 + B w2 H ’ 


”2 * C w2 M 2 > 


x 0 = (0 u f a q x T z T 6T 6th 6s) 
Xs e e 


(30) 

(31) 

(32) 

(33) 
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= (6e 6s 6th) 


(34) 


w 0 = (u* a q ) 
l w w 


X ■ % V X 0 \ X 0 0 0) • <36) 

where A^, B^, and are matrices of appropriate size corresponding to the 
coefficients in the original equations. Expressions for the elements of these 
matrices can be obtained in terms of the stability derivatives of the aircraft. 
Equations (30) - (36) form the design model for the longitudinal control law. 


B. LATERAL DESIGN MODEL 

The design model for the lateral control law is obtained using the same 
approach as the longitudinal (vertical) model. The general equations of motion 
are linearized about a flight condition corresponding to level wings, a constant 
airspeed equal to the desired airspeed, a flight path angle equal to the se- 
lected glideslope angle, with flaps and gear down in the landing configuration. 
The lateral perturbation equations can then be expressed as 


4> = sec0 o (cosY o p 4 siny Q r)_ 


iL = sec0 (sina p 4- eosa r), 
Y o o r o 


6 - a 9 (j) 4 a 99 g 4 a 9/ p 4 a„ r 4 b 91 6A 4 b 99 SR 


4 b 00 Ssp 4 h 01 g 4 h QO p 4 h OQ r 
33 ^ 31 w 32 r w 33 w 


p = a. 9 3 4 a, t p 4 a, _ r 4 b, - SA 4 b. 9 SR + b. 9 Ssp 


4 h 71 6 4 h yo p 4 h, 9 r 

41 w 42 ■ w 43 w 
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( 41 ) 


r = a co 3 + a c/ p 4* a__ r + b_ 1 6A + b co 6R + b co 6sp 

53 54 r 55 51 52 53 r 


+ h 3 + h p + h 

51 w qo 


52 K w 


53 


r 

w 


where 3, P and r are the sideslip angle, roll rate, and yaw rate, respectively, 
3 , p and r are the sideslip angle, the roll rate, and yaw rate due to wind 

WWW 

velocities only, <5A, 6R and 6sp are the perturbations of the ailerons, rudder, 
and spoilers, respectively. The coefficients a^ in the above equations 
depend on the aircraft stability derivatives. 

The position of the aircraft relative to runway centerline is expressed by 
the perpendicular distance of the aircraft center of mass to the runway center- 
line. This distance, normalized by the aircraft r s steady-state airspeed, U q , 
will be used as a state variable in addition to the equations already obtained. 
Then, the lateral distance y (in feet) of the aircraft can be expressed by 

y = U 0 [l es ( 2,1)(1 + u r ) + L es ( 2,2)3 + L ES (2,3)a] , (42) 

where u 1 is the normalized inertial perturbation in the speed along the x 

s 

direction and a is the perturbation in the inertial angle of attack, and 


L E g(2,l) = cosa Q .cos(0 + 8) sinip + cos<|> sina Q sin(0 o + 9) sin^ 

- sina Q sin<j> cos^ (43) 

Lgg (2,2) = sinc() sin(9 Q + 9) sinij; + cos<j) cosip (44) 

L__(2,3) = -sina cos(9 + 9) sinijj + cosa cosd) sin(9 +9) sinil^ 

ES o o r o o 

- cosa o sincj) cosip . (45) 

Rearranging the terms in equation (42) 


( 46 ) 
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where r} 


y 


and y ? 


are given by 


n y = (L es ( 2,2) - l)e + L es (2,1)(1 + U’) - cos(a o - 8 q )* 


+ Lg (2, 3) a + sina Q (j) 


( 47 ) 


y' 


■ yi* o 


(48) 


Note that equation (46) contains no approximation when a in (47) is inter- 
preted as the normalized inertial velocity component in the z g direction, but 
is simply a rearranged form of (42) with the nonlinear terms grouped into a 
single term. To obtain the second order terms for the remaining variables, 
note that 


$ = sec©^ (cosY q p + siny^ r ) + *1^ 

i = sec0 (sina p + cosa r) + n, 
T o o o ^ 

a 

y s 

& = ir ~ r + n 3 » 

where 


(49) 

(50) 


(51) 


n 4> 

n 3 


(cos<b tan(0 + 0) - tan0 ) (cosa r 4- sina 
Y o o o o 

(cosd) sec(0 + 0) - sec0 ) (cosa r + sina 
T o o o o 

pa - ru , 


p) > 

(52) 

p) > 

(53) 


(54) 


when q is assumed to be small. 

The servo systems, hydraulic and mechanical actuator systems on this air- 
craft have clearly noticeable nonlinear effects. These effects include usual 
nonlinearities such as hysteresis, rate and position limits, as well as d.c. 
gain variations as a function of actuator position, dynamic pressure, etc.. 
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and nonlinear spoiler feedforward and feedback systems- These nonlinear 
effects were not included in the !l control design model" described here. How- 
ever, these nonlinearities were included in the simulation model, and design 
changes were made to accommodate these nonlinear effects during the control 
development using the nonlinear simulation and later during the flight tests. 

The controls which affect the lateral motion of the aircraft are the 
aileron, rudder and spoiler surface settings as can be seen from equations (37) 

- (41), where the spoiler action is used mainly to aid the effect of the ailerons 
during turns. Thus, the spoiler setting is programmed according to the aileron 
setting. This is approximated and modeled here as 

Ssp = C <$A , G = 1.73 . (55) 

y spa spa 

A rate command structure is used for the rudder control; hence, the rudder 
position is considered a state variable which is obtained by integrating the 
rudder rate command, 

6R = u 2 , (56) 

where u^ is considered to be the rudder rate control. 

If the relation between the spoiler and aileron given in (55) is substi- 
tuted in (39), (40), and (41), then the spoiler terms are eliminated from the 
equations. Now, forming a state vector x such that 

x T * (<j) ijj g p r y r SR) (57) 

and a control vector u such that 

U T = (iSA SR) , (58) 

equations (37) - (41), (46), (55), and (56) can be combined into a state variable 
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model of the lateral motion of the aicraft of the form 


x = Ax *f Bu 4* Dw + T) , (59) 

T -T 

where w = (3 p r ) and X] = (tj, T). J] 0 0 0 T) 0) . 

www 9 r p y 

The lateral motion of the aircraft is described by the state variable 

model given in equation (59); this model describes the response of the aircraft 

when a control is applied or when the wind velocities such as gusts or steady 

winds are nonzero. The effects of the wind velocities are introduced through 

the vector w. The components of this vector are 3 or the wind velocity along 

w 

the y g direction normalized by the airspeed of the aircraft, p^ or the rotation 

of the air around the aircraft about the x axis, and r or the rotation of the 

s w 

air around the aircraft about the z g axis, respectively. The roll rate p^ and 
yaw rate r^ components of the wind vector w consist only of the effects of wind 
gusts, thus having an average value of zero; i.e. , these components do not have 
a steady state effect but introduce turbulence effects into the equations. On 
the other hand, the $ w or the normalized lateral wind velocity contains terms 
for both wind gusts and steady winds; thus, it is modeled as 


3=3 +3 

w ^g .s 


( 60 ) 


where 3 is the gust or turbulence term, and 3 is the steady wind term. The 
§ s 

gust terms are of a random nature and can be modeled using the well-known 
Dryden spectrum. The transfer functions (i.e., the ratio of the Laplace 
transforms of the gust generating system output to its input) for the gust 
model can be expressed by 


V s) ' -°v 


S i + 


V 


ttU" 


o J 


3 L s 
v 

U 


1 + 


L s' 
v 

U 


(61) 


18 




It should be noted that even though is statistically independent of 

3 w and r , the latter two are not independent of each other. As shown in 

equation (63), the lateral gust, 3 , must be input to the system represented 

§ 

by the transfer function, G (s) , to obtain r . Since r is obtained from 3 , 

J ’ r w w g’ 

the two variables are highly correlated. 

Consider the case of a steady wind. In the Earth- fixed coordinate system, 

the wind velocity has a component in the direction of runway centerline W , 

x 

and a component perpendicular to the runway centerline say W<- ; it is assumed 

that there is no steady wind in the vertical direction although gusts may be 

present. Hence, if L is the transformation matrix from Earth-fixed to sta- 

SE 

bility coordinates, then the steady component, 3 > of the normalized lateral 

s 

wind velocity is given by 

^s = L SE^ 2 ’ 2 ^ W 5 + L SE^ 2,1 ^ W x ’ 
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L_„(2,l) = sin(9 + 9) sin<j> cosib - cos <b sirup 

oHi O 


9 


(65) 


L (2,2) = sin(0 + 9) sin<j> simp + costb cost); . (66) 

oiL O 

To include wind shear into the model, the steady lateral wind velocity 
can be described as 


W 5 = W 6 + ai 3 , (67) 

W 6 = 0) 4 , (68) 


where oo^ and U)^ are Gaussian white noise processes independent of each other 

and of 8 , p and r . Now, the transfer function for the gusts described in 
g r w w & 

(61), (62), and (63) can be combined into a state variable model of fourth order. 
Adding (67) and (68) to this model, we obtain a sixth-order model of the form 

W - A W + B a) , (69) 

w w 

where (3^ can be expressed as 


- B g + 6 S - w i + w 5 + h <70 > 

?1 - (Lgg(2, 2) - 1)W 5 + L se (2,1)W x . (71) 

Defining the remaining elements of the vector £ to be zero, w can be written 
as 


w = C W + 5 (72) 

w 

This expression can now be substituted into (59) to obtain 

x = Ax + Bu + DC W + n (73) 

w 

n = d£ + n (74) 
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The design model for the lateral dynamics, including the aircraft aero- 
dynamics, the actuators and wind conditions, is given by (73) and (74) . 

C. STOCHASTIC SAMPLED-DATA FORMATION 

As the aim of the study is to design and flight test a digital automatic 
landing system obtained by modern control techniques, this section formulates 
an optimal control problem using the design models obtained in the previous 
sections. The solution to this optimal control problem is used as an initial 
design in the development of DIALS. Although some modifications to the design 
were made to account for nonlinearities and other unmodeled effects, the basic 
structure of the optimal control law was left unchanged. 

The optimal control problem is formulated as a stochastic sampled-data 
quadratic regulator with random disturbances [3] - [8] • The longitudinal and 
lateral control laws are treated separately, but using the same basic approach. 
Both the longitudinal and lateral plant models can be expressed in the form 

x — Ax + Bu + DC W + TJ , (75) 

w y 

W = A W 4* B (jo . (76) 

w w 

The trajectory of the aircraft during a final approach and landing includes 
the localizer and glideslope capture and tracking phases with the capability to 
crab and decrab, followed by the flare maneuver. Thus, the model of the desired 
trajectory must be general enough to accommodate the various phases of flight 
involved. It is further desirable to obtain a constant gain control law for the 
various phases of the final approach and landing, to the extent possible. To 
accommodate these characteristics, the desired trajectory is modeled in the form 

z * k z + E , (77) 

z z 
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( 78 ) 



where £ is assumed to be a white noise process. The desired trajectory z(t), 
is thus modeled as a random process whose statistical properties are determined 
by (77), (78). The assumption of a random trajectory is a conceptually appealing 
one, as it stresses the fact that the trajectory may be altered at any time; 
the future values of the trajectory are thus not certain, but can be predicted 
using the model of the trajectory. In general, for an arbitrary deterministic 
trajectory, the optimal control depends on future values of the trajectory [l2], 
except when the trajectory can be expressed as the output of a homogeneous 
dynamic system. 

The error in the actual trajectory, x(t), can now be defined as 

e(t) = x(t) - z(t) (79) 

The models of the aircraft dynamics (75), (76) and the desired trajectory (77), 
(78) can be combined, and expressed in terms of the trajectory error, e(t), 
in the form 


e = Ae + Bu + DC W + (A - A )z - £ 4*p » 

w z * 

W * A W + B 03 , 

w w 

z = A z + r 
z z 

* 

where to, and ^ are assumed to be Gaussian white noise processes. 
Rearranging the terms in (80) - (84), 


(80) 

(81) 

(82) 

(83) 

(84) 
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e = Ae + Bu + Dd , 


(85) 


d = A - d 4- 5 
d 


( 86 ) 


where 


A7 T T r T Tv 
(W z C z T] ) 


(87) 


and ^ is a white noise vector. In (85), (86), the disturbance form of the 
problem is apparent, where the disturbance, d, is uncontrollable by u, and is 
the output of a marginally stable system. 

As the basic objective of the control law is to track the desired trajec- 
tory, a cost function which reduces the error vector, e, is suitable. Consider- 
ations of the closed-loop system damping make it desirable to also include 
cross-coupling terms. As a quadratic objective function accommodates these 
characteristics, results in a linear feedback law and is computationally 
feasible, the objective function is selected as 

J c = ~2t~~ f ^ s|e T (t) Q e(t) + 2 e T (t) M u(t) + u T (t) R u(t)| dt (88) 

To obtain a discrete control law for digital implementation, the usual constraint 
of a constant control between sampling instants is placed 


u ( t ) — u, , kT < t < (k+I)T . (89) 

k o o 

With the assumption of (89), the system equations in (85) and (86) and the 
cost function in (88) can be expressed in terms of the state variables at the 
sampling instants [s] o 

Vi ■ *% + ru k + r d d k + «!k ■ (90> 
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d k+l d k + ? 2k 


(91) 


where e^ and d^ represent the samples e(kT Q ) and d(kT Q ), and are 

white noise sequences. Similarly, the cost function (88) can be discretized 
in the form 


J = 



A T 1 /\ T 1 r P /\ 

Q % + \ R u k + 2 e k N d k + 2 e k M \ 

+ 2 4 s \] 


(92) 


It is important to note that the discrete cost J and the continuous cost 
differ by a constant independent of the control sequence used. Thus, to 
find the control sequence which minimizes the continuous cost J , it suffices 
to obtain the control which minimizes J. The control which minimizes the cost 
J subject to the constraints of (90) and (91), is given by [7], 


u k H lk e k H 2k d k ’ 


(93) 


where e^ and d^ are the least mean-square estimates of e^. and d^, respectively, 
given past and current measurements, and the gains and are given by 


-1 


H lk =, K 
G lk ■ ^ 


Ik 


G lk ’ 

H 2k - G 2k 

9 


(94) 

p lk* + fi ’ 

°2k ■ ^ D k + S 

9 


(95) 

r d + P 2k ^d* 

A A T" 1 

\- E + r p ik 

r , 


(96) 

are given by the nonlinear 

difference equations: 


T T 

■ p ik*-°ik 

A— I A 

*k G ik + G * 

P 1N 

A 

= Q , 

(97) 

* - r ^ G ik 

T D k + « . 

P 2N 

A 

= N , 

(98) 
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Thus, over a finite optimization interval t^ of N sampling intervals, the 
optimal control law is specified by (93) - (98) , It should be noted that the 
gain and the cost matrix P ^ are the optimal solutions to the LQG problem 
without any disturbance; i.e. , d = 0, As the optimization interval N increases, 
it is well-known that -> H^; and P^, under loose conditions [13] , where 

the closed-loop system matrix <j> - FH, is stable. If the disturbance matrix <J)^ 
satisfies p((f)^) < 1, it can be shown that as the optimization interval 

N increases without bounds, and the cost, J, converges to a finite value. It 
should be noted that, due to the random plant noise added at every sample, the 
averaging of the cost over the optimization interval is necessary to maintain 
a finite stochastic cost in any (stochastic) LQG problem. The gains and 
can then be obtained by solving (94) - (98). Using these steady-state gains 
and expanding the components of d^, the control law in (93) can be rewritten as 


U k=- H e e k - H w W k - H Z Z k - H ? S zk - % \ 


It is seen that the control law feeds back the estimates of the error in 
the trajectory to stabilize the error dynamics with the same gain that would be 
obtained from the no wind, no path following case. However, the control law 
uses additive terms depending on the desired trajectory, z^ and 5^, for tracking 
purposes, and uses the second order nonlinear term rj^ as well. Finally, the 
wind estimate W^_ is used directly in the control law to prevent offsets due to 
wind conditions. The estimates of the trajectory error and wind velocities are 
obtained using steady-state Kalman filters driven by the onboard sensors and 
the MLS signal. 

The phases of the final approach and landing considered are the localizer 
and glide-slope capture and track, crab/decrab and flare to touchdown. The range 
of flight conditions considered is 2.5° - 5.5° for glideslope angle, 115 - 135 
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knots for airspeed* 70*000 - 90*000 lb for aircraft weight* .2 - .3 for c.g. 
location* and crosswinds below 15 knots. 

The main considerations in the development of DIALS include the overall 
smoothness of the flight during the final approach, low overshoots of the 
localizer and glideslope, quick settling on the trajectory and accurate tracking 
in adverse weather conditions including gusts and shear winds. However, the 
relative importance of these desirable characteristics varies with the phase of 
the final approach considered. For example, low overshoot characteristics under 
differing wind conditions are particularly important during the localizer and 
glideslope capture maneuvers whereas the actual capture path is of secondary 
importance; however, during the flare maneuver, the vertical path error is more 
significant, and, in particular, the touchdown pitch attitude must be suffi- 
ciently positive to avoid landing on the nose gear. Similarly, when tracking 
the localizer and glideslope, it is important to avoid any offsets and have low 
sensitivity to wind gusts. A more detailed description of the desirable char- 
acteristics in the various control modes is given in Section IV, To achieve 
these varying objectives, it appears that some changes in the control law, from 
one phase of the flight to the next, are necessary. 

Some of the differing objectives for each phase of the final approach can 
be accommodated by different values of the weighting matrices (Q, R) in the 
cost function. However, in an optimal control design, this would necessitate 
changing all of the gains of the control law at the initiation of each mode, 
along with possible transient problems due to switching gains. To avoid pos- 
sible complexity, the approach taken was to obtain a basic design considering 
the total trajectory, and then vary a small number of gains using easy-on f s 
and to modify the structure by switching integrators in or out appropriately, 
at the initiation of each mode. Thus, the basic design model does not contain 
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any integrators; however, when the objective requires it, integrators are 
switched in or out. For example, at the initiation of the localizer (or glide- 
slope) track mode, the integral of the lateral (or vertical) offset is introduced 
to obtain a type 1 system without steady offsets. Thus, the basic design is 
obtained as a stable closed-loop system to which modifications are switched in 
or out at each mode according to the main objectives of that mode. The specific 
modifications made at each phase are detailed in Section IV. 

Nonlinearities in the system and unmodeled effects in the aerodynamic 
characteristics, actuator systems, sensors arid electronics were simulated in 
as much detail as available to obtain a realistic computer simulation of the 
overall system. As these nonlinearities are not included in the design model, 
the optimal control law obtained does not directly accommodate their effects. 
Although the robustness of the control law will indirectly accommodate some 
of the nonlinearities and unmodeled dynamics, the remaining effects on the 
performance were then reduced or eliminated by minor modifications to the con- 
trol law, using the simulation. 

Thus, the approach to the overall design of the digital automatic landing 
system (before flight tests) was: 

1) obtain an optimal control law for a nominal condition as a starting 
point, 

2) using realistic nonlinear simulations, modify and update this law to 
account for the various nonlinearities in the actuator servomechanisms, hydraulic 
systems, engine dynamics and the aerodynamic response, as well as to enhance 
desirable characteristics which may not adequately be reflected in a quadratic 
cost function, 

3) identify a small number of modifications which can be made at each 
phase of flight to achieve the differing objectives of each mode. 
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It should also be noted that the cost function ( 88 ) penalizes the error 
in the actual trajectory relative to the desired trajectory; however, it also 
penalizes the total control u. As noted in [ 7 ] , such a cost function does not 
necessarily provide sufficient trim unless the cost matrices are appropriately 
selected or the control cost is modified to penalize the deviation from a pre- 
selected trim [14]. For example, consider the action of the elevator and stabi- 
lizer. On the TSRV B-737 research aircraft, the effects of the elevator and 
stabilizer on the motion of the aircraft are the same, except that the effec- 
tiveness of the larger stabilizer surface is twice that of the elevator. It 
is desirable to make trimming changes using the stabilizer, andf maintain the 
elevator near its equilibrium position (referenced at 0 ) to avoid high moments 
on the elevator hinge over sustained periods of time. Of course, the fast 
motion (due to the actuators) of the elevator is used to stabilize the dynamic 
modes of the aircraft. Since the trim values required on different glideslopes, 
or wind conditions, or c.g. locations, etc. vary with the particular condition 
encountered, it is desirable to have the stabilizer position move to trim the 
aircraft, while the elevator is maintained at low hinge moment values except 
in transients. Thus, it is desirable not to penalize stabilizer position, but 
to penalize elevator position to obtain the desired trimming action. Since, 
it is also desirable (and necessary) to have slow stabilizer motion, high 
stabilizer rate should be appropriately penalized to obtain smoother transi- 
tions. Thus, appropriate selection of the cost matrices with rate command 
structure can provide appropriate trimming with optimal control formulation 
described in the previous section. Similarly, throttle position is not penal- 
ized, while a nonzero penalty on throttle rate provides smoother and slower 
autothrottle action. In general, control activity can be placed at desirable 
levels by penalizing the control rate. 
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To obtain a 3-D control law, the variable corresponding to distance along 
runway centerline (x^) was eliminated from the equations of motion for the 
design model. However, the filter development model included (and estimated) 
this variable. 

D. STATE AND JWIND ESTIMATION 

The sampled-data formulation of the design problem requires estimates 
of all the state variables in its implementation. As the states corresponding 
to the desired trajectory are known, these variables are not estimated. However, 
the aircraft motion variables and, in particular, the wind velocity variables 
including steady winds, shear winds and gusts are required for quick control 
adjustments to changing wind conditions in the full-state feedback formulation*. 

The MLS position information and the on-board sensor signals are fed to 
the longitudinal and lateral DIALS filters shown in Figures 3 and 4. The fil- 
tered estimates are then used to compute the control surface commands as well 
as to define the desired path. The form of filter used is a discrete Kalman 
filter with constant gains. The steady-state Kalman filter gains were obtained 
as a starting point, then some gains were modified to obtain better performance 
using a nonlinear simulation of the aircraft. The use of constant gains for 
the Kalman filters reduces the considerable computation requirements imposed 
by the error covariance updates. The use of discrete filters along with digital 
controls is well suited to the discrete nature of MLS guidance system. The goals 
of the filter development include the accurate estimation of the aircraft’s 
position, velocity, attitude and wind velocities using on-board sensors usually 
available on commercial aircraft, but without using costly inertial platforms, 


*The author’s recent development of a fast, reliable and convergent algorithm 
for the stochastic output (limited state) feedback problem has made the use 
of optimal dynamic compensation along with complementary filters practical, 
and more desirable in high order problems [15] . 
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and angle of attack or sideslip sensors which are not currently available on 
many aircraft,, 

The aircraft’s position is obtained using the MLS guidance system, which 
provides volumetric coverage in the terminal area. The aircraft receives range, 
azimuth and elevation information at discrete intervals from which it can obtain 
its position with high accuracy even under adverse weather conditions. The 
ground azimuth antenna is located at the runway centerline with coverage up to 
±60°. The DME antenna which provides the range of the aircraft is generally 
co- located with the azimuth antenna. If the DME is located to the side of the 
runway, a simple transformation can be used to obtain the aircraft’s position. 
The elevation antenna is located at the glidepath intercept point (GPIP), but 
is offset to the side of the runway; it provides the aircraft’s elevation angle 
up to 20°. Thus, the aircraft has accurate position information in the volume 
of space within the limits mentioned above. 

Consider a right-handed coordinate frame with its origin at the phase 
center of the azimuth antenna, the x-axis along runway centerline and positive 
towards the runway, and the z-axis positive vertically upwards. If the posi- 
tion of the aircraft in this coordinate system (the MLS coordinate frame) is 
(x o , y Q , z q ) , then the MLS signals have the values given by the formulas below. 

R = \nc 4- y^ + z^ (100) 

* o J o o 


Az = sin 


-1 



( 101 ) 


El = tan — . : ===== srsa== =s- (102) 

V (x o - x E )2 + (y Q - y E > 2 

where x^ and y^ are the x and y coordinates of the elevation antenna phase 
center in the MLS coordinate frame. 
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The onboard sensors used are three body-mounted acceleromaters , attitude 
gyros for the pitch, roll and yaw angles, attitude rate gyros, baromatric al- 
titude and sink rate, airspeed and radar altitude. The filter does not require 
expensive inertial platforms for sensor measurements. All of the onboard sen- 
sors are usually available on commercial aircraft except for the body acceler- 
ometers. However, the acceleromaters are relatively inexpensive instruments and 
in many of the newer aircraft, a normal accelerometer already exists for pitch 
axis control. It should be noted that the radar altitude measurement is used 
in place of the MLS elevation signal only after the aircraft crosses the runway 
threshold since the aircraft flies out of the MLS elevation coverage during the 
flare maneuver. 

To provide accurate estimates of the wind velocity components in the 
longitudinal, lateral and vertical directions, the filters were formulated using 
the aerodynamic properties of the aircraft. In addition, estimates of the bias 
errors in the three body-mounted accelerometers, barometric altitude and sink 
rate and attitude sensors are also obtained by the filters. The lateral filter 
uses the roll and yaw attitude, yaw rate and roll rate, the processed MLS y^ 
measurement, and the lateral accelerometer. The remaining measurements are 
used for the longitudinal filter. Discrete longitudinal and lateral Kalman 
filters were formulated with steady-state gains, with some cross-coupling using 
second order terms and wind estimates. The measurement bias terms were not 
included in the Kalman filter formulation and the steady-state gain computations, 
but the required gain matrices were obtained using simple approximations. The 
basic form of the filters can be expressed by 
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(105) 



(106) 

(107) 

(108) 


where x^, W^, b^ are predicted values and x^, W^, r|^ are filtered values. 

It should be noted that whereas the control gains are based on the longi- 
tudinal and lateral sampled-data formulation of the previous section (i.e., 

(90) and (91)), the design of the Kalman filters do not use the same models in 
obtaining the filter gains. As most of the variables in d^ are known trajectory 
parameters and a (nonlinear) function of x^, only attitude, attitude rate, 

position, inertial velocity and wind velocity need to be estimated or simply be 
filtered for noise. Thus, in the above equations, the control position variables 
in the vector x^ are not included in the estimation, and is computed using x^. 

The position equations used in the filter implementation have a strong 
dependence on the glideslope angle, y , selected for the particular approach 
from the continuous equations (9), (10). While the filter gains F , F and F, 
are kept constant, the elements of the <J> matrix corresponding to the position 
variables are modified to account for different glideslope angles. Since the 
glideslope angle has to be entered into the flight computers prior to glideslope 
capture, the computation of these elements can be done in real-time. 

Table 1 shows the sensor error statistics used in a nonlinear simulation 
of the sensor outputs. The filter designs were tested in simulation prior to 
the DIALS flight tests. During the simulation, some gains were modified to ob- 
tain better performance of the overall control law. 
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III. IMPLEMENTATION OF THE DIGITAL CONTROL LAW 


The digital integrated automatic landing systems (DIALS) is designed for 
flight testing on the Transport Systems Research Vehicle (TSRV) , a Boeing 
737-100. The research aircraft is equipped with various onboard sensors, four 
General Electric 703 digital computers*, as well as monitoring, testing, and 
recording facilities. The digital control law is implemented on the GE 703 
flight computers. 

The flight computers are 16 bit, fractional machines with fixed-point 
arithmetic. The constraint of fixed-point arithmetic with 16 bit words makes 
the software implementation of a complex digital control law considerably more 
complicated than a floating-point arithmetic machine. To insure the accuracy 
of the control command computations, it is necessary to meticulously keep track 
of every variable in the computer memory for range and resolution. Since most 
variables in the software are intermediate variables, not explicitly treated in 
the design of the control law, this task results in considerable effort. 

A. SCALING AND ROUND-OFF ERROR ANALYSIS 

By considering the physical meaning of the variables involved, and the 
expected flight profile, safe upper and lower limits for each variable can be 
assigned. These limits determine the scaling necessary for each variable. As 
long as the resolution with 16 bits does not violate the accuracy requirements 
for that variable, a single precision representation may be possible. Otherwise, 
a double precision representation is necessary for that variable. As the flight 
computers use two’s complement arithmetic an overflow can have devastating ef- 
fects, and must be guarded against with great caution. From the point of view 

*The GE 703 computers are currently being replaced by Norden 1170 computers. 
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of overflow* only the position variables require double precision representation 
in the implementation of DIALS . 

The accuracy of the control law computations also depends on round-off 
errors introduced in the multiply-add operations required. The computations 
required to obtain the control commands are described by (99), the filter com- 
putations by (103) “ (108). Equation (99) results in the commanded control 
variables, and is obtained by a finite number of multiply-add operations whose 
accuracy is easy to analyze. However, the filter computations are recursive, 
so that errors introduced in previous computations are processed through the 
filter every iteration, and may produce much larger errors. If the errors 
introduced due to round-off exceed the accuracy desired, it is necessary to use 
double precision arithmetic to reduce these errors. 

An analysis of the accumulation of the round-off errors can be made in the 
following manner. The filter equations in (103) - (108) can be expressed in 
the form 


= Yi " C X. , 
k •'k ok 3 


\ \ + F o \ ’ 

\+l - *0 \ + F o u k ’ 


(109) 

( 110 ) 
( 111 ) 


where and correspond to predicted and filtered variables, respectively, 
and are stored in the same locations in computer memory. Each of the variables 
in the vectors X^., U^, and y^ are scaled by the diagonal matrices $ x , S^, 
and S , respectively. 

Let the superscript 11 ^ " denote the actual representation of a vector in 
the computer. Then the operations actually performed can be expressed in the 
form 
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(112) 


(113) 


(114) 


where e^, and are the round-off errors introduced in performing the 
operations implied by (112), (113), and (114), respectively. These errors 
depend on the order in which the multiply-add operations are performed, the 
number of bits retained in a product when the add operation is performed, and 
the actual number of operations which produce round-off error. In general, 
these errors are nonlinear funtions of the variables which are complicated to 
treat analytically. However, under some assumptions [16], [17] , the statistical 
nature of the error can be modeled as a sequence of uncorrelated random vectors 
with uniform distribution. When a large number of operations are involved the 
probability distribution function resembles the Gaussian distribution [l 7] . 

Denoting the scaled coefficient matrices in (112) - (114) by the subscript 
"l", these equations can be combined to obtain 


K+l ■ [+1 - *1 F 1 C l] *k + (+1 F 1 F l) F k + r ! C k + 


(115) 


£, = £a, + <}>,£, + d), F, e . 

k xk r l xk r l 1 Vk 


(116) 


The representation X^ corresponds to in engineering units, thus, let 


E k = \ - S X \ > 


(117) 


[ <J> - <j> F C E, +£, , 

|_ r o r o o oj k k 


(118) 


e, = <j> F S £ , + T S e , + S £, 
k o o y yk o u uk x k 


(119) 
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The matrices n^* n x ^> n V k are diagonal, with the entry in each row corresponding 
to the number of times a round-off error is introduced in the corresponding 
equation in (112) - (114). Usually this is the number of multiply operations 
performed in each scalar equation. 

It should be noted that the effects of scaling, the order in which the 
equations are implemented (i.e., (112) - (114)), the word size, as well as some 
special procedures can be computed and analyzed with this approach. Tables 2-7 
shows the resolution, range and accuracy required for the HALS implementation. 

When appropriate scaling is done, the error introduced by a small number 
of round-off errors can usually be neglected. From (99), it can be seen that 
the error in the control commands u^ due to round-off can be approximated by 
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/s •J* 

H E t with covariance H Ra H , 
e k e J3 e 

Consideration of these effects shows that the position computations, and 
the innovation computations should be performed in double precision to have 
sufficient accuracy in the control commands. 

B. CONTROL IMPLEMENTATION FORM 

The Experimental Avionics System Integration Laboratory (EASILY) at the 
Langley Research Center (LRC) provides a facility where a control law implementa- 
tion can be tested in real-time. The basic procedure is to code the control law 
to be tested in the flight computer which will be used onboard the aircraft. The 
flight computer is then brought to the EASILY lab and interfaced with a detailed, 
nonlinear, six-degrees-of-f reedom, real-time digital simulation of the aircraft 
and its subsystems. The facility contains CRT displays which simulate the dis- 
plays on the TSRV research aircraft so that the simulated flight can be monitored. 
As the digital real-time simulation of the aircraft dynamics and its subsystems 
is located in a different building, the flight computer and displays at the 
EASILY laboratory are connected to the simulation through communication lines. 

While the communication lines linking the digital aircraft simulation to 
the flight computer containing the control law are expected to introduce some 
delay in the signals, the magnitude of such delays was not originally estab- 
lished. Tests of the end-to-end delays show that these can reach 200 msec, 
with an average lag in the vicinity of 100 msec. The large magnitude of this 
delay is partially due to the signal conditioning at both the EASILY lab and 
at the digital simulation facility. A further difficulty is introduced by the 
fact that the digital flight control computer and the digital aircraft simulation 
computer operate asynchronously. This introduces delays of random magnitude, 
partially explaining the variation in the measured amounts. 
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Considering the fact that DIALS operates at a low sampling rate of 10 Hz (a 
sampling interval of 100 msec*) ? delays of the same magnitude clearly introduce 
significant difficulties in its performance* It is well-known that delays as 
low as half the sampling interval of a control law can destabilize the closed- 
loop system [lg.]* Thus, the actual implementation of DIALS was modified to 
accommodate delays of the magnitude of 100 msec* 

To minimize the software and design impact, the following approach is used 
to implement the DIALS equations. Consider the system equations given in (90) 
and (91), and the form of the control law given in (93). Let p £ 0, and consider 

u k - p Vi ■ - H i( 5 k - p Vi) " H 2(v p 5 k-i) • (124) 


From the filter equations, it can be seen that 


e k * e k-i + r d d k-i + r u k-i + G i \ 


d k " d k-l + G 2 A 


(125) 

(126) 


Substituting (125) and (126) into (124) , regrouping terms neglecting 
results in 

"k - ( p 1 - H i r )vr “il 4 ’ - p "Ki- [ H 2 (+d- p r a] d k-i (127) 

It is important to note that in this form, u^ depends only on previous 
samples, and can be computed ahead of time; so that if the simulation has a 
known delay, then u^ can be commanded as soon as it is computed to counteract, 
or reduce, the actual delay. While during the flight tests, it can be commanded 
at? its intended time. It is important to note that neglecting the innovation 
terms corresponds precisely to using the predicted state rather than the filtered 
state; i.e.. 
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(128) 


u 


k 




Also note that, except for the first term, (127) is of the same form as 

before in terms of software, whereas e, is obtained at a later time, Further- 

k 

more, (127) introduces a smoothing of the control commands which is preferable. 

It is interesting to note that the closed-loop poles with this control 
form are the same as before. However, new poles have been introduced, as the 
order of the system has increased by augmenting the state by u^. The new system 
contains a real closed-loop pole at A = p of order r, the number of controls. 
This can be seen by noting that the eigenvalue/ eigenvector equations 


({> e + F u = X e , 

(129) 

- - p l)e + (p I - H r r)u = X u , 

(130) 

hold if, and only if, 


(j)e + ru = Xe, andju = - K e or A = . 

(131) 

Now, if u equals - K e, then 


U - r k) e = X e , 

(132) 


so that X is an eigenvalue of ^ <}> - F K^, the closed-loop system with (93). 
While if X = p, p is an eigenvalue of the new system. The eigenvalues equal 
to p have linearly independent eigenvectors. Appropriate choice of p, i.e., 
-1 < p < 1, maintains the stability of the closed-loop system. 

C. MISCELLANEOUS IMPLEMENTATION ISSUES 

While it is almost impossible to mention all the issues that come up in 

designing and implementing a complex system such as DIALS, a few of these are 
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mentioned in the following. 

1. Stabilizer Trim 

DIALS is designed with an independent stabilizer trim control. The 
stabilizer command has the tendency to seek a stabilizer position where the 
elevator position is zero. This is desirable to reduce the moment load on 
the elevator hinges. 

On the other hand, the stabilizer trim design on the TSRV is difficult 
to modify to allow independent operation. Thus, the implementation of DIALS 
in the flight tests uses the TSRV’s original trim to reduce costly mechanical 
modifications to the aircraft systems. 

The approach used to modify the original DIALS design is to form a new 
elevator command by combining the original elevator and stabilizer commands. 
Noting that the stabilizer effectiveness in the flight regime considered is 
approximately twice as high as that of the elevator, the elevator command 
implemented adds the DIALS stabilizer command, multiplied by a factor of two, 
to the original elevator command. This elevator command and the TSRV stabilizer 
trim function operate simultaneously. 

The stability of the system with this modification was investigated along 
heuristic lines, by noting that the elevator command implemented produces 
essentially the same moments on the aircraft as the original design. The sta- 
bilizer trim only acts to rearrange the proportion of these moments produced 
by the elevator and the stabilizer. 

2. Sampling Rate 

The design of the control and filter gains in DIALS is based on a sampling 
rate of 10 Hz (i.e., 100 msec sampling interval) for both the control commands 
and the sensor signals. On the other hand, the GE 703 digital flight computers 
on which DIALS is implemented has a major cycle of 49.152 msec. Thus, DIALS 
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is implemented at a rate of 2 major cycles per iteration, or a sampling interval 
of about 98,3 msec, rather than the 100 msec used in the design. 

To assess the magnitude of the effects introduced by such a mismatch in 
the sampling rate, consider 

<p(T+e) = e = e e ^ e (I+Ae) , (133) 

r(T+e) = / o T+e <Kt) dx B = T(T) + (J) (T) r(e) * T(T) + <f) (T) Be . (134) 

It can be seen that the effect of a sampling mismatch can be viewed in 
terms of controlling a modified plant with gains designed for a somewhat dif- 
ferent plant. The variation in the plant parameters is of first order in the 
sampling interval mismatch £, with sensitivity directions along <J>(T.).A and <J>(T)B. 

During the development phase of DIALS, in numerous simulations corresponding 
to various flight conditions, the original gains have been slightly modified 
to attain a level of insensitivity to parameter variations. Thus, it would be 
expected that, given the small magnitude of the mismatch, the effects would not 
be large. Simulations of this condition show that the effects are not partic- 
ularly noticeable when sensor noise and gusts are present. The flight test 
results tend to confirm the. simulation results. 

3. Use of "Easy-on 11 Smoothing 

During the final approach and landing, the aircraft captures the localizer 
and glideslope (sometimes simultaneously), tracks the localizer and glideslope, 
crabs and decrabs, and flares to touchdown. Each of these phases introduces a 
new mode in the control law, as will be described in the following sections. 

In each mode, the commanded trajectory changes. Similarly, when the aircraft 
reaches a low altitude, the MLS elevation signal degrades and is essentially 
unavailable past its antenna located at the glidepath intercept point (GPIP). 
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Thus 3 the DIALS filter switches from the elevation signal to the radar altimeter 
signal. 

To enhance the smoothness of the flight and avoid undesirable transient 
effects, easy-on smoothing is used in switching from one commanded path to 
another, when sensor signals are switched, or when control gains are changed. 

The easy-ons used are mostly linear and exponential. The period over which an 
easy-on operates is obtained by heuristic considerations, simulation results 
and flight tests, as the theoretical formulation does not include this particular 
type of smoothing. The effect of the smoothing is to avoid transients excited 
by sudden changes in velocity and acceleration; the stability of the system, 
however, is unaffected by such smoothing. 

4. Barometric Sink Rate 

The longitudinal Kalman filter in DIALS uses a sink rate measurement 
obtained from barometric altitude augmented by a normal accelerometer. During 
the implementation phase, it was found that the barometric sink rate signal 
available on the aircraft is obtained by mechanically differentiating the 
barometric altitude signal, and does not contain a normal accelerometer. So 
that this sink rate signal is of very low quality, particularly in turbulence, 
with errors reaching 30 - 40 ft /sec depending on the wind conditions. 

The implementation of DIALS during the flight tests uses a sink rate 
signal obtained by complementary filtering the MLS altitude signal with body- 
mounted accelerometers avoiding a redesign of the longitudinal Kalman filters. 

5. Engine Pressure Ratio Versus Thrust 

While the thrust provided by the engines at a given level of throttle is 
not directly measured, the engine pressure ratio (EPR) for both engines is 
measured and available in the flight computers. An estimate of the thrust can 
be obtained using the measured values of the engine pressure ratios. 
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To obtain an analytical relationship between EPR and thrust for the range 
of expected flight conditions, data are obtained using steady state or trim 
values of EPR and thrust from a detailed simulation of the aircraft and its 
subsystems for a range of mach numbers varying from .17 to .24. A c.g. location 
of .25, altitude of 1000 ft, and various glidepath angles, airspeeds, weights, 
etc., are used in obtaining the data. A linear fit of the form 


T = a (EPR) + b 


(135) 

a = 26.7878, b = - 25.9395 

y 

(136) 


fits the data obtained with a maximum error of 30 lb over the range of mach 
numbers expected. This relationship is used to obtain an estimate of the thrust 
is the units of 1000 lb. As the data are obtained from the simulation model, 
and for steady-state conditions, the actual accuracy of the estimate in a 
dynamic, gusty environment is not clearly determined. 

Finally, the flight software was tested successfully in static tests, by 
comparing results with those obtained from the simulation of DIALS on the LRC 
Computer System. Dynamic tests of the flight software were also performed. 
Dynamic tests require the use of an aircraft simulation as part of the dynamics. 
Due to synchronization requirements, the aircraft simulation used in the dy- 
namic tests was implemented on a different computer than the one used for con- 
trol law design. Thus, due to differences in the two aircraft simulations, 
the dynamic tests were only able to produce a less-than-perf ect check of the 
flight software. Despite lengthy debugging attempts, some software bugs made 
their way to the flight tests and were discovered during the flights, while 
others probably were not. 
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IV. FLIGHT TEST RESULTS 


The Digital Integrated Automatic Landing System (DIALS) described in the 
previous sections was flight tested on the Transport Systems Research Vehicle 
(TSRV) , a Boeing 737-100, to test and demonstrate a modern control direct- 
digital-design automatic landing system. The flight tests were conducted by 
the National Aeronautics and Space Administration’s Langley Research Center at 
Wallops Island Center. This section discusses the performance of the automatic 
control system during each phase of the final approach and landing and describes 
the corresponding control mode in further detail. 

The phases of flight for which the flight tests were conducted follow. 

1) Localizer capture 

2) Steep glideslope capture 

3) Localizer track 

4) Glideslope track 

5) Decrab 

6) Flare 

The localizer and glideslope captures can be performed independently or 
simultaneously depending on the aircraft position relative to the Glide Path 
Intercept Point (GPIP). Another feature of the automatic landing system is 
that it can perform steep approaches. This is achieved by selecting the glide- 
slope angle at any time prior to glideslope capture within the limits of 2.5° 
to 5.5°. Typically, the aircraft is positioned, manually or automatically, be- 
low the selected glideslope so as to intercept the localizer at a selected air- 
speed, as shoVn in Figure 1. When the localizer or glideslope capture criteria 
are satisfied, the corresponding capture mode is engaged; so that depending on 
the initial aircraft position and attitude, the localizer and glideslope can 
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be captured simultaneously or sequentially. As the aircraft reaches the local- 
izer or glideslope, the corresponding track (or hold) modes are engaged. When 
a cross-wind component is present, the control system crabs the aircraft into 
the wind, followed by a decrab maneuver when a decrab altitude is reached. A 
flare path which depends on the glides lope selected for a given approach is 
generated on-line, and the aircraft is controlled about this path until touch- 
down. 

The runway at NASA Wallops Flight Center used for the DIALS flight tests 
is equipped with a high accuracy Microwave Landing System (MLS). The MLS re- 
ceiver onboard demodulates the signals and feeds the elevation azimuth and DME 
information into the flight computer. 

A. LOCALIZER CAPTURE 

In the localizer capture mode, the control law tries to perform a smooth 
transition from the aircraft’s initial position to the runway centerline, so 
that the aircraft can quickly settle and stabilize on the localizer. Two 
characteristics are stressed in the design and development of DIALS for this 
mode: the reduction of overshoots, and rolling away from the localizer during 

the capture. Consistently low overshoots of the localizer are particularly 
important in airports with close parallel runways. This property can be used 
to increase the efficiency of the terminal area as it allows the two parallel 
runways to be operated independently. A low overshoot reduces the settling 
time on. the localizer, thus allowing close-in captures when necessary. 

The second charaeteristis is to initiate the capture maneuver by rolling 
the aircraft away from the runway centerline. Pilots usually initiate localizer 
capture in this fashion, and consider rolling towards the runway centerline 
"the wrong direction". This property can be obtained by initiating the local- 
izer capture when the natural tendency of the control law is to roll in the 
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desired direction. To achieve this property, consider the component of u, in 

K 

(99) corresponding to the aileron command. The distance, y , from the local- 

K 

izer at which the aileron command crosses zero was selected to initiate the 
localizer capture mode. Assuming level wings and using the inital conditions 
for the filter states, and neglecting small terms, a localizer capture criterion 
of the following form can be obtained 


* U ol c l \ + °2 \ + c 3 B wk 
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where y^ 3^, and 3^ are estimates of the distance from runway centerline, 

the heading relative to the runway, the inertial normalized side velocity and 
the normalized cross-wind velocity, respectively. 

The localizer capture mode is initiated automatically when the inequality 
holds. The constants are chosen so that the initial tendency is to roll away 
from the localizer with the possible exception of conditions when the aircraft 
is already rolled as desired, say due to gusting winds. It should be noted 
that the localizer capture mode is engaged farther away from the runway center- 
line when the desired speed is higher, as would be expected. To further 
ensure a smooth transition, an "easy-on" function was used on the aileron 
command signal. 

To avoid large overshoots beyond the runway centerline, appropriate 
closed-loop damping of the control law was obtained by proper selection of 
the weighting matrices Q and R, with subsequent testing in the simulation of 
the localizer capture. 
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Figures 5 - 10 show the localizer capture phase of the landing during 
the flight tests* As can be seen from the plots, the localizer overshoot is 
rather small even in the presence wind conditions with 15 knot cross-wind com- 
ponents. The flights test included localizer intercept angles up to 50°. The 
average value of the localizer overshoot over all the flights is 24.2 ft, with 
a standard deviation of 25.7 ft, as shown in Table 8. 

It should be noted that the aircraft can be considered to have settled 
on the localizer within approximately 40 sec in most cases. Also note that 
the initial tendency of the aircraft is to roll away from the runway, as desired. 
It should also be noted that the design does not use an integrator on the side- 
slip angle in order to null it; so that during the capture the sideslip angle 
is not centered about zero. Although in steady wind conditions the angle is 
small, in high gust conditions it can reach higher values. It is desirable to 
modify the design to produce lower sideslip, particularly in high gust condi- 
tions, without modifying the other characteristics of the control law. 

It should be noted that while all three runs show the presence of cross- 
winds, the flight shown in Figure 9 corresponds to high turbulence conditions, 
as can be seen by the cross-wind estimate and the variation in the attitude. 

B . GL IDESLOPE CAPTURE 

The desirable properties of the glideslope capture mode include low over- 
shoot and the tendency to pitch down at the initiation of the capture. To have 
a low overshoot of the glideslope for the various steep glideslope angles that 
may be selected (i.e., 2.5° - 5.5°), it is desirable to initiate the capture 
when the aircraft is well below the glideslope. This approach is different 
than the practice of engaging the glideslope hold mode when the glideslope is 
intercepted (e.g., see[9j), where an overshoot, even if small, cannot be 
avoided. 
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Initiation of the capture mode below the glideslope alone does not guarantee 
a low overshoot. Thus, some experimentation in the selection of the weighting 
matrices coupled with simulation was used to achieve this property. In partic- 
ular, for both overshoot characteristics and smoothness of maneuver, the sink 
rate command was eased on from zero to the glideslope sink rate. This imple- 
mentation provides a simple method of adjusting the speed with which the capture 
maneuver is performed without an appreciable overshoot. 

While engaging the capture mode below the glideslope is desirable for 
overshoot characteristics, it also produces an initial tendency for the aircraft 
to pitch up. To avoid this effect, the capture engage logic was selected so 
that the initial tendency of the elevator is to produce a negative pitching 
moment . 

To achieve this property, consider the component of u^, in (99), corres- 
ponding to the elevator command. If the glideslope capture mode is engaged at 
the altitude, h^g^, where the elevator command crosses zero, it can be seen 
that the initial elevator command will tend to gradually produce a negative 
pitching moment. 

While the thrust also has a considerable contribution to the pitching 
moment in this aircraft, due to the engine location, the DIALS autothrottle 
and engine dynamics have much longer time constants than those associated with 
the elevator. Thus, changes in the throttle command and thrust are of small 
magnitude over a short period of time, so that the initial pitching tendency 
is mainly determined by the faster elevator command. Although the initial 
tendency of the DIALS throttle command is to lower the thrust level, which is 
the desirable action, the glideslope capture criterion was selected by consid- 
ering only the elevator action. After some manipulation, and neglecting terms 
of small magnitude, the following glideslope capture altitude can be obtained. 
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(140) 

(141) 

(142) 


where Y Q is the selected glides lope angle, x^ is the estimated ground distance 

to the glide path intercept point (GPIP) , e^ is the estimated error in the state 

x^, and is the commanded inertial angle of attack. Note, that, as in the 

localizer capture, the glideslope capture altitude is lower when the desired 

speed is higher. Further note that the criterion is applicable to the capture 

of an arbitrary glideslope by using the glideslope angle selected by the pilot 

as y in (140) . 
o 

Figures 11 - 16 show the glideslope capture phase of the final approach 
during the flight tests. The glideslope angle desired for the final approach 
is entered into the flight computer prior to the glideslope capture. Until the 
glideslope capture criterion is met, the longitudunal dynamics are controlled 
usually by a different control system designed using classical methods. 

When the glideslope capture criterion is met, the system automatically 
engages the glideslope capture mode. In a given approach the point at which 
capture is initiated varies largely according to the initial position of the 
aircraft, the glideslope selected, and the airspeed selected. 

On the plots shown in Figures 11, 13, and 15, it is seen that at glideslope 
capture initiation, the control law produces the desired combined actions of 
pitching the aircraft down even though it is still below the glideslope, and 


49 



lowering the thrust level* Accordingly, the sink rate is brought to the level 
necessary to track the glideslope selected for the approach under the pre- 
valent wind conditions, at the selected airspeed. 

It should be noted that depending on wind conditions, and the performance 
of the control system in operation prior to the capture, the aircraft pitch 
rate, sink rate, and thrust rate may be noticeably different than zero, even 
though a steady level flight is desired. The flight results indicate an 
insensitivity to such variations, presumably due to the capture criterion. 

So that the capture is initiated in a coordinated fashion under the various 
wind conditions encountered. 

The approach shown in Figure 13 corresponds to a conventional 3° glide- 
slope, whereas Figures 11 and 15 show steep glideslope captures corresponding 
to a 4.5° glideslope. The flight results indicate that the control law per- 
formance is insensitive to the particular glideslope selected. Except for 
the pitch attitude and sink rate levels needed for the particular glideslope, 
the basic behavior of the longitudinal variables during the capture appears 
to be similar, variations being largely due to wind conditions. 

As seen in the plots of vertical path deviation, the glideslope over- 
shoots are quite small for steep as well as conventional glideslopes. The 
average glideslope overshoot in the flights performed is 4.5 ft, with a 
standard deviation of 2.3 ft, as shown in Table 9. In fact, it may be noted 
that the overshoot is no larger than the vertical path deviations on the 
glideslope which are due to wind gust disturbance. 

The time required to settle on the glideslope selected is seen to be 
about 20 sec from a performance standpoint, quick settling on the glideslope 
is desirable as it enhances close-in captures. On the other hand, from a 
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passenger comfort standpoint , a slower , more gradual capture may be pre- 
ferable. 

Finally , it may be noted that the glideslope capture maneuver appears to 
produce no undesired transient activity in the airspeed, as measured by the 
calibrated airspeed (CAS). So that, even through the aircraft is pitching 
down during capture, the reduction of the thrust during the same period 
largely counteracts the effects of the pitch on airspeed while enhancing the 
pitch-down action. The latter effect is due to the fact that the engines are 
located below the aircraft center-of-gravity , thus, reducing thrust also re- 
duces the pitching moment. 

C. LOCALIZER AND GLIDESLOPE TRACK 

As the aircraft gets closer to the localizer or glideslope, the local- 
izer or the glideslope track (or hold) modes are engaged, the order depending 
on the particular approach path. In comparison to the capture modes, the 
main objectives considered in the design of the track modes were to achieve 
quick settling on the localizer and glideslope, insensitivity to wind dis- 
turbances (gusts, shear and steady winds), and accurate tracking of the 
localizer and glideslope. 

Since it is not desirable to start the flare maneuver before the air- 
craft is fully stabilized on the glideslope, a low settling time for the 
track modes is important to enable "close-in captures". The capability to 
accurately capture and track a shallow or steep glideslope provides a flex- 
ibility which can be used to avoid vortices when following large aircraft, 
reduce the noise perceived by airport communities, and generally increase 
terminal area operations efficiency. 
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The localizer track mode is engaged when the aircraft is 30 ft from 
the runway centerline or 25 sec after the localizer capture has been engaged. 
Since the main objective of the track mode is to minimize excursions away 
from the runway centerline, the gains feeding back the cross-track error esti- 
mate to the aileron and rudder are increased gradually for a smooth mode 
transition. Since a type 1 system is desirable in this mode, an integrator 
for the cross-track error is also initiated to avoid steady -state offsets 
from the runway centerline. Initially, the integrator gains were set equal 
to a tenth of the gains for the cross-track error, i.e. , the product of 
the sampling interval and the cross- track error /gain. The gains were then 
adjusted using simulation results to improve the performance under the 
various conditions expected to occur. 

The glideslope track mode is engaged when the aircraft intercepts the 
glideslope , or 25 sec after the glideslope capture mode has been engaged. 

The latter criterion is used to ensure that the track mode is engaged even 
if the aircraft does not overshoot and cross the glideslope. To reduce any 
excursions from the selected glideslope, the gains feeding back the deviation 
from the glideslope and the error in the normalized vertical velocity are 
increased gradually using an easy-on. To avoid steady-state offsets from 
the glideslope, an integrator on the glideslope deviation is introduced 
when the track mode is initiated. 

During both the capture and track modes, the DIALS filters estimate 
the components of the wind velocities. These estimates are fed to the con- 
troller to compute the surface command signals. Thus, to the extent that 
the wind velocities are estimated, the control law performs a gust allevia- 
tion function; i.e., even if no position or inertial velocity errors are 
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present, the wind velocity estimates will produce surface commands which tend 
to counteract the effects of winds, thus reducing excursions away from the 
glideslope. The wind velocity estimates are also used to maintain a con- 
stant airspeed. This is obtained by commanding an inertial speed deviation 
equal to the negative of the wind estimate in the along- track direction. The 
commanded airspeed thus remains constant at the desired nominal value. 
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The performance of the system in the localizer and glideslope track 
modes are shown in Figures 5-16. As mentioned in the previous sections, 
the aircraft settling time on both the localizer and glideslope is short. Due 
to the lack of overshoots, by the time the track modes engage, the aircraft 
is largely settled on both the localizer and glideslope. 

The main function of tracking the localizer and glideslope is per- 
formed satisfactorily, with small deviations about the desired path. In the 
high gust case, the localizer offsets during the track mode are seen to be 
below 30 ft while the aircraft adjusts to gusts and changing wind condi- 
tions. The variations in the roll and yaw angles also appear to be appro- 
priate for the high gust conditions encountered. 

The glideslope deviations in the high gust condition tend to be of the 
order of 8 ft except for a rather high deviation of about 16 ft presumably 
due to a high vertical gust. In general, the sink rate and pitch behavior 
seems appropriate. However, the thrust is seen to have an oscillatory motion, 
which may be due to unmodeled throttle or engine dynamics. 
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Figure 16(a) shows the throttle rate command of the control law and the 
measured throttle position* Comparing the zero crossings of the throttle 
rate with maxima and minima of the measured throttle position shows that 
the throttle lags its command by 1 - 2 see. As this large lag, presumably 
in the throttle servo, is not included in the design model the control law 
does not account for its existence in the control loop. As the nonlinear 
simulation used in the development of the control law does not contain 
such a long time constant in the throttle servo, the effect was not noticed 
in the simulations, and only became clear in the flight tests under high 
gust conditions where the throttle rate commands are high and may reach 
their rate limit. 

D. A WIND SHEAR ENCOUNTER 

An excellent example of the performance of the control law under wind 
shear conditions was encountered on the flight shown in Figure 7. As seen 
from the yaw angle behavior, the aircraft, has a crab angle while tracking 
the localizer due to an estimated cross-wind velocity of about 12 knots. 

At about 100 sec, the crab angle is seen to be reduced to almost zero. At 
first glance, this appears as a normal decrab maneuver. However, the decrab 
altitude has not been reached. 

The reason for the behavior can be seen in the estimated cross-wind 
velocity, which indicates that the cross-wind of approximately 12 knots has 
very quickly died down, and may even have reversed direction. Thus, the 
behavior of the yaw angle was a response to a sudden change in the steady 
wind. Similar adjustments are seen in the roll angle and localizer deviation. 
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Although a wind change of about 15 knots in 15 sec (or approximately 
150 ft of altitude on a 3° glideslope) which corresponds to a wind shear of 
about 9 knots per 100 ft is not as severe as some reported cases, the ease 
with which DIALS handles this particular encounter demonstrates a desirable 
performance. The magnitude of the maximum localizer deviation is seen to 
be less than 30 ft. The flight continued to make a normal automatic landing. 

E, CRAB AND DECRAB 

In the presents of cross-winds, the desired behavior is to crab into 
the wind, i.e., to head into the wind with level wings and no slidslip while 
remaining on the localizer. If this condition is not specifically accommo- 
dated, an optimal quadratic regulator will usually both yaw and roll into 
the wind, thus having a nonzero steady-state sideslip angle which is not de- 
sirable. In addition, the bank angle can be uncomfortable. To obtain the 
crab condition in DIALS, a roll integrator is fed back when the aircraft roll 
reaches 2°, provided that the localizer track mode has been engaged. This 
ensures a zero steady-state roll angle forcing the yaw angle to automatically 
adjust itself to the angle necessary for crab. To obtain crab with a quick 
and smooth transient response, the commanded (aerodynamic) sideslip angle is 
set to zero , i.e. , 



(145) 


where z^ is the normalized lateral velocity (inertial sideslip) and B ws 
is the estimate of the steady wind velocity. As a small sideslip is desir- 
able in all the phases of flight, except for decrab, the commanded sideslip 
angle is set to zero at localizer capture initiation and remains active 
until decrab. 
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During decrab the aircraft heading is aligned with the runway center- 
line 9 while the roll angle is used to maintain the aircraft on the localizer. 
The decrab mode is initiated at an estimated altitude of 250 ft. It should 
be noted that this decrab altitude is somewhat arbitrarily selected, and 
can be reduced to initiate the maneuver at a lower altitude. To obtain 
the decrab maneuver, first the commands which produce the crab condition 
are phased out; the roll integrator is phased out and the inertial side- 
slip command z^ is set to zero gradually. Since this corresponds to a 
nonzero commanded aerodynamic sideslip angle, the wind velocity estimates 
fed into the controller produce a decrab tendency. To enhance this ten- 
dency and ensure decrab as a steady-state condition, an integrator of the 
heading relative to the runway is introduced to aid the aileron and rudder 
actions necessary for decrab. Initially the integrator gains were set 
equal to a tenth of the heading error gains. These were later adjusted using 
simulation results to improve the performance throughout the flight regime. 

Figures 5 - 10 show the aircraft lateral dynamics during crab and 
decrab. After the localizer capture, the roll angle settles at a average 
value of zero as desired until the decrab altitude- Similarly, the sideslip 
angle remains small. On the other hand, the heading relative to the runway 
settles at a nonzero value in order to maintain the localizer deviation at 
zero. It may be noted that the yaw angle settles at the crab angle quickly 
and smoothly. This is essentially due to the use of cross-wind velocity 
estimates in the control law. To the extent that the estimated cross-wind 
is correct, the yaw angle settles smoothly at its desired value, and the 
localizer overshoot is small. As the control law does not use a measure- 
ment of the sideslip angle, this information is provided by the wind velocity 
estimate. 
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When the decrab altitude is reached, the decrab mode is engaged. During 
the flight shown in Figure 5, the cross-wind estimate shows a steadily de- 
creasing cross-wind velocity, along with a steadily decreasing crab angle. 

So that when the decrab mode is engaged, the aircraft is nearly aligned with 
the runway. 

On the other hand, the flight shown in Figure 9, the cross-wind has 
a steady component of about 20 knots throughout the final approach, gusting 
to over 20 knots on occasion. The decrab mode reduces the yaw angle to 
align the aircraft heading with the runway. The decrab law contains a roll 
limit of 4.5°; i.e. , the system decrabs as much as possible up to a 4.5° 
bank angle, as rolling further entails greater risk of the engine touching 
the runway if a heavy gust were to occur at touchdown. Thus, any remaining 
misalignment would be taken by. the landing gear, or the landing would be 
aborted by the pilot in higher cross-winds. As can be seen, the performance 
of both crab and decrab appears to be satisfactory in these flights. 

F. FLARE 

During the flare maneuver, the main objectives are to reduce the air- 
craft’s sink rate to an appropriately lower level and to touch down with a 
sufficient margin on the pitch attitude, aligned with the runway, near the 
selected touch-down point. As DIALS is required to flare from steeper glide- 
slopes than usual for a B-737, the aircraft has a higher sink rate than usual 
on the higher glideslopes. The approach taken was to generate a flare tra- 
jectory on line as a function of the glideslope angle, the desired touch- 
down flight path angle and the touch-down point; so that when the glide- 
slope is steeper, the flare (initiation) altitude is higher. It should be 
noted that, for a given approach, the flare path is fixed in space, and is 
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fed to the automatic control system as the desired or commanded trajectory . 

The family of flare paths was developed starting from the desired 
vertical acceleration profile. The following form was selected for this 
purpose 
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AX f = (P - 2 X f ) 12 , (148) 

where is the x-coordinate of the aircraft’s position vector, X^ the 
touchdown point, the altitude at X^ (zero for touchdown), Y q and y^ 
the selected flight path angles on the glides lope and at touchdown, respec- 
tively. The desirable characteristics of this profile included the smooth- 
ness of transition (e.g., h ff , hence h, has a continuous first derivative at 
flare initiation), the simplicity of the parametric form (i.e. , flare profiles 
from various glideslopes are obtained according to the value of y ) , and 
the ease with which h" can be integrated analytically to obtain an expression 
for h f and h to determine an altitude profile for flare. 

It is of interest to note that the flare profile described above was 
selected after some experimentation with an altitude versus x^ profile of 


*The prime 11 1 " denotes differentiation with respect to the distance 
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exponential form. These formulations resulted in undesirable transient 
behavior at flare initiation, presumably due to the discontinuity in the 
commanded vertical acceleration when starting flare. The cosine type ac- 
celeration profile provides a smooth transition in this variable, while 
incorporating all the parameters of interest in the flare path such as 
glideslope angle, touchdown point, etc. For a constant ground speed, x^, 
the vertical acceleration profile is seen to be 


h (x ) = x h * (x ) 
e e e 


(149) 


h(x ) = x 2 h"(x ) (150) 

e e e 

• •• 

The vertical profile h(x^), Mx^), Mx^) is used in the desired path or 
command vectors, z^ and in order to follow the flare path for the se- 
lected glideslope. However, the altitude profile alone does not guarantee 
all of the desirable and critical properties necessary for an acceptable 
touchdown. The pitch attitude at touchdown is of utmost importance. A 
negative pitch angle at touchdown would result in the nose gear touching 
down before the landing gear and being subjected to high levels of loading. 
Thus, a positive pitch angle safety margin is usually desirable. It is 
also usual to have a decreasing airspeed profile during flare to aid the 
pitch profile during flare, as well as ease the touchdown speed, sink rate,, 
and rollout. To obtain these characteristics the command variables and 
were set as follows. 
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where ^ is the desired flight path angle at time t^, GEZ5 is an easy-on, 

Z* , is the desired altitude in the Earth coordinate axes, L. , is the 
ed k xj 

estimate of L (i,j), and e 1 is the first column of the identity matrix. 

Eb 1 

DELVFR and G are constants which can be set according to the amount of 
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decrease in airspeed and the pitch profile desired, respectively. These 
values were also modified using the profiles obtained in simulation runs. 

The performance of the automatic landing system during the flare maneuver 
is shown in Figures 11 - 16 . It should be noted that the plots of the vertical 
path deviation in these figures show the difference between the aircraft al- 
titude estimate and the flare path generated on-line corresponding to the 
selected glideslope. The on-line flare path was generated with the equations 
shown above. 

During the DIALS flights, ten ?, hands-off M automatic landings were per- 
formed. Of these, three approaches were made with a selected glideslope of 
3°, while steep approaches on 4.5° glideslopes were selected for seven flights. 

The plots of the sink rate during flare show that the control law reduces 
the sink rate from its level on the glideslope to a lower level for touch- 
down. This is largely achieved by pitching up to increase the lift, as well 
as to obtain a positive pitch attitude at touchdown. The vertical path de- 
viation remains small throughout flare; however, the deviations at touchdown 
largely determine the actual touchdown point. 

Of the ten automatic landings made, the flight where the wind shear was 
encountered resulted in the highest touchdown sink rate of 3.5 ft /sec, which 
is quite acceptable even if considered hard. The other two landings resulted 
in touchdown sink rates of 2 ft/sec, and 3 ft/sec. The average touchdown sink 
rate of all the automatic landings was 2.4 ft/sec, with a standard deviation 
of 0.74 ft/sec. 

It should be noted that the commanded touchdown sink rate at 125 knots 
is 2.2 ft/sec. Thus, the touchdown sink rate statistics with its low standard 
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deviation shows the desired performance. The average touchdown point was 
1026 ft, with a standard deviation of 244 ft. The commanded touchdown point 
was 1152 ft. Tables 10 and 11 show the touchdown statistics for the automatic 
landings performed. 

The control performance during flare shows similar characteristics for 3° 
and 4.5° glideslopes. The flare altitude for a 3° glideslope was 102 ft, 
with an altitude of 134 ft for a 4.5° glideslope. While the flare altitudes 
used are higher than usual, a higher than usual flare altitude is reasonable 
for steeper glideslopes. Nevertheless, the particular values used do not 
have a special significance, and can be easily modified by an appropriate 
modification of flare path parameter values. 

G. FILTER PERFORMANCE 

In the context of control law design, the filter plays the role of a 
dynamic compensator; i.e. , a dynamic system whose outputs are available 
for use in controlling the plant. In the sampled data LQG problem posed 
in section II C, as well as other LQG formulations, the optimal dynamic 
compensator is the least-mean square estimator of all the states. Thus, 
the filter smooths noisy measurements and produces estimates of the vari- 
ables which are not directly measured. The smoothed signal necessarily 
lags the measured signal. The use of the dynamic plant model in the filter 
provides important lead information which is used to reduce the lag. However, 
if the accuracy of the model used is not sufficient, significant lags can 
be produced which tend to destabilize the overall closed-loop system, and 
may lead to the divergence of the filter itself. 
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As discussed earlier* DIALS uses fixed-gain (steady-state) Kalman 
filters for both the longitudinal and lateral dynamics. As the flight 
conditions considered vary from one approach to the next, the plant dynamics 
model used in the design is certain to have mismatches in some flight con- 
dition. Thus, it is important to accommodate errors in the filter design, 
model appropriately. This entails the use of plant noise to account for 
model uncertainties, and has the effect of increasing the filter gains at 
the expense of some extra noise in the estimates. 

A feature of DIALS is the direct use of wind velocity estimates in 
the control law. In particular, the estimation of cross-wind velocity is 
necessary for use in the crab and decrab maneuvers, but more significantly, 
in reducing the localizer overshoot. While the inertial cross- track velocity 
can be estimated using the MLS position and body-mounted accelerometers, a 
measurement of the lateral airspeed, such as a sideslip sensor, was not al- 
lowed in the design, although a sideslip sensor is available on the TSRV. 

So that the main sensor available for wind estimation is the airspeed sensor, 
where the cross-wind is of second order. Thus, to estimate the cross-wind 
speed with some accuracy, and more significantly, to estimate sudden changes 
in the cross-wind without long lags, it is necessary to use more information. 
The aircraft dynamic model provides such information even though the model 
may have some parameter uncertainty. 

From the flight performance of the overall automatic control system 
discussed in the previous sections, it is clear that, the filters are pro- 
viding the necessary dynamic compensation without introducing large lags in 
the feedback loop, despite the presence of plant model errors. Comparison 
of the longitudinal wind estimate with the airspeed signal in Figures 11, 
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13, and 15 indicates a small lag which may be improved. However, it is 
also clear that the filter, appropriately, does not attribute every 
variation in the airspeed to a change in wind speed, and does not show the 
initial tendency to do so and correct itself later. 

As the actual lateral wind velocity is not measured, the accuracy of 
the cross-wind estimate cannot be directly ascertained. However, the 
"sure-footed" manner of the crab decrab maneuvers while remaining on the 
localizer indicate that the error in the estimate must be small. More 
importantly, the response to the wind shear condition shown in Figure 7 
indicates that the lags in the estimate are small, although perceptible. 
Comparison with the sideslip angle is also interesting. 

Plots of the innovations for Flight No. 2 which contains the wind shear 
condition discussed are shown in Figure 17, It is interesting to note that 
the noise in lateral position and altitude get smaller as the aircraft 
approaches the runway, and therefore the MLS antennas. As these variables 
are obtained from the angular measurements of azimuth and elevation, the 
noise variance is directly proportional to the range of the aircraft. Note 
that the altitude innovations show a transient behavior which stabilizes 
about zero. Also note that the roll and pitch innovations show a small bias. 

In general, the innovations show a behavior indicative of white noise 
which is the expected behavior. However, considering the known modeling 
uncertainties, some coloring of the innovations should be expected. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The flight testing of the Digital Integrated Automatic Landing System 
(DIALS) is among the currently rare occasions where a complex modern control 
digital design has been demonstrated in flight. To the author's knowledge, 
it is the first automatic landing system designed by optimal direct-digital- 
design techniques which has been successfully flight tested. As indicated 
in the report, the design of a modern control system for a complex, flight 
critical task is not straight-forward, despite a large body of literature on 
some aspects of the techniques involved. In the following, some areas of 
theoretical research which address problems encountered in designing practi- 
cal automatic control systems are also recommended. 

DIALS has six modes corresponding to the phases of the final approach 
and landing: (1) localizer capture, (2) steep glidelsope capture, (3) local- 

izer track, (4) glideslope track, (5) crab/decrab, and (6) flare. The control 
law uses MLS position information, body-mounted acceleromaters, and sensors 
usually available on commercial jets, but does not require inertial platforms. 
Using the MLS, DIALS has the capability of performing steep approaches to 
landing by selecting the desired glideslope angle. 

During the DIALS flight tests, ten "hands-of f 11 automatic landings were 
performed on the TSRV aircraft. These automatic landings were performed 
under diverse wind conditions, including a wind shear case. The final ap- 
proach flights included the modes of localizer capture, steep glideslope 
capture, localizer and glideslope track, crab/decrab, and flare to touchdown. 
On three of the approaches, a conventional 3° glideslope was selected. On 
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the remaining seven automatic landings, steep approaches with a 4.5° 
glideslope were made. 

All the modes of the control law were exercised during the flight tests. 

The flight test results indicate improved performance in the following areas 
under all wind conditions encountered: (1) reduced localizer overshoots, 

(2) reduced glideslope overshoots, (3) reduced settling time for both localizer 
and glideslope, (4) flexibility of performing steep final approaches, (5) accu- 
rate tracking of the localizer and glideslope, (6) smooth crab and decrab maneu- 
vers, (7) smooth touchdowns from steep glideslopes, and (8) excellent handling 
of a wind shear condition encountered. It should be noted that the accuracy 
of the localizer and glideslope tracking relative to the runway is enhanced by 
the high accuracy of the MLS signals. It should also be noted that the sampling 
rate of the control system is 10 Hz, which is lower than that usually required 
in analog designs which are implemented digitally. In general, the DIALS 
flight tests demonstrate that modern control methods can successfully be used 
to design automatic control laws for complex tasks. 

Areas where the performance of the control law can be improved include: 

(1) the inclusion of the lags present in the throttle loop into the design 
model to improve auto throttle performance and (2) the reduction to the side- 
slip angle during localizer capture. The importance of the design model is 
clearly illustrated by the autothrottle performance. Of course, classical 
design techniques also require adequate modeling of the plant, as evidenced 
by similar autothrottle behavior observed in flight tests of a classical 
control system on the same aircraft. As a simulation, no matter how detailed, 
is only an approximation of the actual plant, the importance of flight testing 
is also illustrated. 
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DIALS is a full-state feedback law as it is based on a sampled-data 
LQG formulation. For numerous reasons, including flexibility of control 
structure and ease of implementation, it is desirable to design control 
laws where the variables used in the feedback loop can be selected by the 
designer. Thus, the development and application of output feedback and 
dynamic compensation provide useful avenues of further research. * 

The insensitivity of a control law to particular parameter variations, 
as well as possible variations in all the parameters has always played an 
important role in control system design. In the development of DIALS, in- 
sensitivity to flight condition parameters such as glideslope angle, air- 
speed, weight, c.g. location, steady wind conditions, sensor jitter, 
sampling rate variations, etc., was obtained by extensive simulations in 
a heuristic manner. However, this is a lengthy and difficult approach. 
Research on design methods where insensitivity, or robustness, are a part 
of the cost function can provide solutions which shorten the design time 
and effort considerably. 

Finally, practical control laws are designed for systems which most 
often contain nonlinearities. Few systems are free of nonlinearities such 
as hysteresis or backlash, nonlinear gains, etc. Such well-known nonlin- 
earities often produce undesired instabilities in control designs which 
do not accommodate them. It would be desirable to develop modern control 
methods and control structures which accommodate such often encountered 
system nonlinearities. 
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FIGURE 1. DIALS FLIGHT PATH GEOMETRY 
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FIGURE 3. DIALS FUNCTIONAL BLOCK DIAGRAM 



FIGURE 4. BLOCK DIAGRAM OF FEEDBACK LOOP 
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FIGURE 7(a). LATERAL FLIGHT PATH VARIABLES (FLIGHT NO. l) 
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FIGURE 7(c). LATERAL FLIGHT PATH VARIABLES (FLIGHT NO. l) 
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FIGURE 9(b). LATERAL FLIGHT PATH VARIABLES (FLIGHT NO. 3) 
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FIGURE 9(c). LATERAL FLIGHT PATH VARIABLES (FLIGHT NO. 3) 
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FIGURE 11(a) . LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. l) 








FIGURE 11(b). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. l) 
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FIGURE 11(c). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. l) 
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FIGURE 12(a). LONGITUDINAL CONTROL VARIABLES (FLIGHT NO. l) 
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FIGURE 12(b). LONGITUDINAL CONTROL VARIABLES (FLIGHT NO. l) 
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FIGURE 13(a). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. 2 ) 







FIGURE 13(b). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. 2) 
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FIGURE 13(c). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. 2) 
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FIGURE 14(b). LONGITUDINAL CONTROL VARIABLES (FLIGHT NO. 2 ) 



FIGURE 15(a). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. 3.) 
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FIGURE 15(b). LONGITUDINAL FLIGHT PATH VARIABLES (FLIGHT NO. 3) 
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FIGURE 16(a). LONGITUDINAL CONTROL VARIABLES (FLIGHT NO. 3) 
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FIGURE 16(b). LONGITUDINAL CONTROL VARIABLES (FLIGHT NO. 3) 
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FIGURE 17(a). POSITION INNOVATIONS (FLIGHT NO. 
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FIGURE 17(b). ATTITUDE INNOVATIONS (FLIGHT NO. 2). 
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TABLE 1. SENSOR ERROR MODEL PARAMETERS 


Sensor 

Noise standard 
deviation 

Bias 

Comments 





Attitude gyro 

.229° 

.229° 


Rate gyro 

.02° /sec 

0 


Body-mounted 

accelerometers 

1% of g 

1% of g 

.229° misalignment 
and 

.25% scaling errors 
also modeled 

MLS - aximuth 

.01° 

.0125° 

jitter added 

MLS - elevation 

j 

.01° ' 

.0125° 

jitter adddd 

MLS - DME 

2.29m 

2. 29m 


Radar altimeter 

. 305m 

. 305m 


Barometric altimeter 

. 305m 

1.52m 


Barometric sink rate 

. 305m/ sec 

. 6 lm/ sec 


Airspeed indicator 

2% 

0 

Noise is multiplicative 


2 

g = 9.81 m/sec 


102 








L03 














TABLE 3. RANGE/UNITS FOR MEASUREMENT VECTORS 


Variable 

Min. 

Max. 

Units 





y Ll 

-0.5 

0.5 

* rad 

y L2 

-0.5 

0.5 

rad/ see 

y L3 

-550.0 

100.0 

sec 

y L4 

-55.0 

1.0 

see 

y L5 

- 55.0 

1.0 

see 

y L6 

-.5 

0.5 

non- dim. 

y L7 

-.125 

.125 

-1 

sec 

V! 

t* 4 

00 

-.5 

.5 

non-dim. 

y L9 

-.125 

.125 

-1 

sec 

y l 

-1.0 

1.0 

rad 

y 2 

-IT 

7T 

rad 

y 3 

-0.5 

0.5 

rad/ sec 

y 4 

-0.5 

0.5 

rad /see 

y 5 

-32. 

32. 

sec 

y 6 

-.125 

.125 

-1 

sec 
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TABLE 4, RANGE/UNITS AND 



FOR ESTIMATED AND PREDICTED STATE VECTORS 


Units 

Computational Accuracy 



rad 

8.7 x 10~ 5 rad (.005 deg) 

non- dim. 

2.3 x 10 4 (.05 ft/sec) 

rad 

5.2 x 10~ 4 rad (.03 deg) 

rad /sec 

6.1 x 10 3 rad/sec (.0035 deg/sec) 

sec 

1.6 x 10 -4 sec (.035 ft) 

sec 

2.7 x 10 -4 sec (.06 ft) 

1000 lbs. 

5 x 10" 3 /1000 lbs. (5 lbs) 

degrees 

8.4 x 10" 3 deg (5 lbs) 

rad 

1.2 x lO" 4 rad (.007 deg) 

rad 

8.7 x 10 4 rad (.05 deg) 

rad 

8.7 x 10" 4 rad (.05 deg) 

rad 

5 x 10 4 rad (.029 deg) 

rad/sec 

N. A. 

rad/sec 

N. A. 

sec 

6.2 x 10 -4 sec (.16 ft) 

rad 

N. A, 



TABLE 5. RANGE/UNITS AND ACCURACY FOR WIND STATES AND BIASES 


Variable 

Min. 

Max. 

Units 

Computational Accuracy 






W LI, W L1 

-.25 

.25 

non-dim. 

2,4 x 10 ^ (.06 ft/sec) 

W L2, W L2 

-.25 

.25 

rad/ sec 

8.8 x 10 5 rad/sec (.022 ft/sec^) 

W L3, W L3 

-.25 

.25 

rad /see 

2.6 x 10 ^ rad/sec (.015 deg/ sec) 

W L4, W L4 

-.5 

.5 

non- dim. 

4 x 10 ^ (.10 ft/ sec) 

W L5, W L5 

-.5 

.5 

non-dim. 

4 x 10 4 (.10 ft /sec) 

W L6, W L6 

-.25 

.26 

rad 

2 x 10 4 rad (.05 ft /sec) 

W L7, W L7 

-.25 

.25 

-i 

sec 

/l o 

1.2 x 10 see (.03 ft/sec ) 

\l j 

-.03125 

.03125 

rad 

1.2 x 10 -4 rad (.007 deg) 

S\ 

b L5 

-2.0 

2.0 

sec 

2 x 10 4 sec (.05 ft) 

A 

b L6 

-.03125 

.03125 

non-dim. 

1.2 x 10 -4 (.029 ft/sec) 

A 

b L7 

-.015625 

.015625 

: 

-1 

sec 

3.2 x 10 -5 sec" 1 (.008 ft/sec 2 ) 

A 

b L9 

-.015625 

.015625 

: 

sec * j 

2.8 x 10 -5 sec" 1 (.007 ft/sec 2 ) 

"i, "1 

1 -.5 

.5 

non-dim. 

-4 

8. x 10 (.2 ft/sec) 

"2, ”2 

-.25 

.25 

-1 

sec 

8. x 10 b sec * (.02 ft/sec 2 ) 

W 3, W 3 

-.25 

.25 

rad /sec 

■ 

1.7 x 10 4 rad/sec (.01 deg/sec) 

A 

w. w. 
4, 4 

-.25 

.25 

rad/sec 

1.7 x 10 4 rad/sec (.01 deg/sec) 

“5, "5 

-.5 

.5 

non- dim. 

8 x 10 4 (.2 ft/sec) 

A 

W, W £ 

6, 6 

-.25 

.25 

non-dim. 

8. x 10 (.02 ft/sec) 

b. 

4> 

b * 

V 

-.03125 

.03125 

rad 

3.6 x 10~ 4 (.02 deg) 

-.03125 

.03125 

rad 

3.5 x 10~ 4 (.02 deg) 

-.015625 

.015625 

-1 

sec 

2.4 x 10~ 5 (.006 ft/sec 2 ) 
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TABLE 7. RANGE/UNITS AND ACCURACY FOR COMMAND VARIABLES 


Variable 

Min. 

Max. 

Units 

Computational Accuracy 






Z L1 

-.25 

.25 

rad 

8.7 x 10 rad (.00r deg) 

Z L2 

-.5 

.5 

non- dim. 

2.3 x 10 4 (.05 ft/sec) 

Z L3 

-.25 

.25 

rad 

5.2 x 10 4 rad (.03 deg) 

Z L4 

-.5 

.5 

rad/sec 

6.1 x 10 ^ rad/sec (.0035 deg/sec) 

Z L6 

-1.0 

1.0 

sec 

1.6 x 10 4 sec (.035 ft) 

? L2 

-.25 

.25 

non-dim. 

2.3 x 10 4 (.05 ft/sec) 

S L4 

-.125 

.125 

rad/ sec 

6.1 x 10 ^ rad/sec (.0035 deg/sec) 

? L6 

-.25 

.25 

sec 

1.6 x 10 4 sec (.035 ft) 

Z 1 

-1.0 

1.0 

rad 

8.7 x 10~ 5 rad (.005 deg) 

Z 2 

-7T 

TT 

rad 

8.7 x 10 ^ rad (.005 deg) 

Z 3 

-tt/4 

tt/4 

rad 

5 x 10 -4 rad (.029 deg) 

Z 4 

-.5 

.5 

rad/sec 

6.1 x 10 ^ rad/sec (.0035 deg/sec) 

Z 5 

-.5 

.5 

i 

rad/sec 

6.1 x 10 rad/sec (.0035 deg/sec) 

Z 6 

-32.0 

.1 

32.0 

sec 

6.2 x 10 4 sec (.16 ft) 

? ul 

-.125 

.125 

rad 

6.1 x 10“ 5 rad (.0035 deg) 

C u2 

-.125 

.125 

rad 

6.1 x 10“ 5 rad (.0035 deg) 

S u3 

-.5 

.5 

rad 

6.1 x 10 -5 rad (.0035 deg) 

^u4 

-1.0 

1.0 

rad/ sec 

_5 

6.1 x 10 rad/sec (.0035 deg/sec) 

C u5 

-.5 

.5 

rad/sec 

6.1 x 10 rad/sec (.0035 deg/sec) 

^u6 

-.5 

.5 

sec 

2.7 x 10 4 sec (.06 ft) 
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TABLE 8. LOCALIZER OVERSHOOT STATISTICS 


Localizer Overshoot (ft) 

Average 

Standard Deviation 



24.2 

25.7 


TABLE 9. GLIDESLOPE OVERSHOOT STATISTICS 


Glideslope Overshoot (ft) 

Average 

S t and ar d Devia t ion 

4.6 

2.3 


TABLE 10. AUTOMATIC TOUCHDOWN STANDARD DEVIATIONS 


Sink Rate (ft /sec ) 

- 1 

Touchdown Point (ft ) 

Average 

Standard Deviation 


Average 

Standard Deviation 






2.4 

0.74 


1026 

244 
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TABLE 11. AUTOMATIC TOUCHDOWN STATISTICS 



Commanded Sink Rate at 125 knots: 2.21 ft/sec. 
Commanded Touchdown Point: 1152.5 ft. 


m 
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